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ABSTRACT 


Computer  programs  are  given  for  calculating  the  characteristic 
cut  rents  and  characteristic  gain  patterns  of  conducting  bodies  of 
revolution.  Also  given  are  computer  programs  for  using  these  char¬ 
acteristic  currents  in  aperture  radiation  and  plane-v/ave  scattering 
problems.  Plot  programs  for  use  with  a  Calcomp  plotter  are  includec ♦ 
Operating  procedures  and  program  details  are  discussed,  and  sample 
input-output  data  are  given. 
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I,  INTRODUCTION 

The  general  theory  and  method  of  computation  of  characteristic  modes 
for  conducting  uodies  of  arbitrary  shape  is  given  in  reference  [ij-  The 
notation  of  this  report  is  consistent  with  that  of  [1],  which  should  be  re¬ 
ferred  to  tor  further  identification  of  the  symbols  used.  The  programs  given 
here  are  those  used  for  the  numerical  results  presented  in  reference  [1]. 

Six  computer  programs  have  been  written.  They  are  defined  according 
to  their  function: 

1.  Calculate  the  generalized  impedance  matrix  Z 

2.  Calculate  the  characteristic  currents  (eigencurrents), 

3.  Plot  the  eigencurrents. 

4.  Calculate  and  plot  the  gain  patterns  of  the  eigencurrents 

5.  Calculate  and  plot  c / *  (scattering  cross  section  divided  by  the 
square  of  the  wavelength)  for  an  axially  incident  plane  wave 

6.  Calculate  and  plot  the  gain  pattern  for  radiation  from  an 
axially  symmetric  excitation. 

These  programs  are  discussed  and  listed  in  the  next  six  sections  Operating 
instructions  and  sample  input-output  data  are  also  given. 

II.  GENERALIZED  IMPEDANCE  MATRIX 

Program  //I  calculates  the  elements  of  the  generalized  impedance  matrix 
for  a  body  of  revolution,  and  is  a  slight  modification  of  the  computer  program 
appearing  in  Appendix  A  of  reference  [1].  The  subroutine  LINEQ  and  its  call 
statement  81  have  been  removed  in  order  to  obtain  the  generalized  impedance 
matrix  Z  instead  of  the  generalized  admittance  matrix  Except  for  NPKI  which 
was  taken  to  be  20  in  all  of  our  work,  all  the  punched  card  data  required  by 
program  #1  is  explained  in  reference  [2],  As  on  page  26  of  reference  12  j, 
punched  ' erd  data  is  cead  early  in  the  main  program  according  to 

50  PJAD(l,5i,END  =  52)  NN,  NP,  NPHI,  BK 

51  FORMAT (3i3,  Eia.7) 

READ (1,5 3) (RH(I) 4  1*1,  NP) 

READ (1,5 3) (ZH(I) ,  I  -  1,  NP) 

FORMAT (10F8- 4) 


53 


Here,  BK  is  the  propagation  constant  k  =  wuie*.  An  odd  number  NP  of  data 

points  are  taken  from  the  generating  curve  C  of  the  body  of  revolution. 

RH(I)  and  ZH(I)  are  respectively  the  distance  p  from  the  axis  of  the  body 

t  h 

of  revolution  and  the  corresponding  z  coordinate  at  the  I  data  point, 

t*  h 

The  first  and  NP  data  points  are  on  the  ends  of  C.  If  C  closes  upon 
itself,  care  must  be  taken  to  make  the  coordinates  of  the  first  data  point 
equal  to  those  of  the  NP^1.  Program  //I  writes  the  NM2  by  NM2  impedance 
matrix 
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on  record  1  of  direct  access  data  set  6.  NM2  is  either  NP-1  or  NP-3  depend¬ 
ing  on  whether  or  not  C  closes  upon  itself.  If  NN  =  0,  zeros  are  supplied 
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I,  INTRODUCTION 

The  general  theory  and  method  of  computation  of  characteristic  modes 
for  conducting  bodies  of  aroitrary  shape  is  given  in  reference  IlJ>  The 
notation  of  this  report  is  consistent  with  that  of  [1],  which  should  be  re¬ 
ferred  to  for  further  identification  of  the  symbols  used.  The  programs  given 
here  are  those  used  for  the  numeri.al  results  presented  m  reference  [1]. 

Six  computer  programs  have  been  written.  They  are  defined  according 
to  their  function: 

1.  Calculate  the  generalized  impedance  matrix  Z 

2*  Calculate  the  characteristic  currents  (eigencurrents) , 

3.  Plot  the  eigencurrents. 

4.  Calculate  and  plot  the  gain  patterns  of  the  eigencurrents - 
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5.  Calculate  and  plot  zh  (scattering  cross  section  divided  by  the 
square  of  the  wavelength)  for  an  axially  incident  plane  wave^ 

6.  Calculate  and  plot  the  gain  pattern  for  radiation  from  an 
axially  symmetric  excitation. 

These  programs  are  discussed  and  listed  in  the  next  six  sections  Operating 
instructions  and  sample  input-output  data  are  also  given. 

II.  GENERALIZED  IMPEDANCE  MATRIX 

Program  //I  calculates  the  elements  of  the  generalized  impedance  matrix 
for  a  body  of  revolution ,  and  is  a  slight  modification  of  the  computer  program 
appearing  in  Appendix  A  of  reference  [2].  The  subroutine  LINEQ  and  its  call 
statement  81  have  been  removed  in  order  to  obtain  the  generalized  impedance 
matrix  Z  instead  of  the  generalized  admittance  matrix  Except  for  NFHI  which 
was  taken  to  be  20  in  all  of  our  work,  all  the  punched  card  data  required  by 
program  //I  is  explained  in  reference  [2].  As  on  page  26  of  reference  [2]* 
punched  card  data  xs  read  early  in  the  main  program  according  to 

30  READ(1 ,5i  ,£ND  =  52)  NN,  NP,  NPHI,  BK 
51  FORMAT (313,  E14.7) 

READ (1,53) (RH(I) ,  1=1,  NP) 

READ (1 ,53) (ZH(I) ,  1=1,  NP) 

53  FORMAT (10F8, 4) 


Iter ef  bK  is  Che  propagation  constant  k  =  w/pe.  An  odd  number  NP  of  data 
points  are  taken  from  the  generating  curve  C  of  the  body  of  revolution* 

RH( 1 )  and  ZH{I)  are  respectively  the  distance  p  from  the  axis  of  the  body 
ot  revolution  and  the  corresponding  z  coordinate  at  the  I data  point* 

The  first  and  NP*'*1  data  points  are  on  the  ends  of  C.  If  C  closes  upon 
itself,  care  must  be  taken  to  make  the  coordinates  of  the  first  data  point 
equal  to  those  of  the  NF^*  Program  //I  writes  the  NM2  by  NM2  impedance 
matrix 

Z  =  Z^  where  n  «  NN  (1) 

n  n  n 

n  n 

on  record  1  of  direct  access  data  set  6.  NM2  is  either  NP-1  or  NP-3  depend¬ 
ing  on  whether  or  not  C  closes  upon  itself.  If  NN  -  0,  zeros  are  supplied 

for  Z^t  and  Z^ . 
o  o 
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Listing  of  Program  01 

II  ( 0034 ,EE  *4, 2 ) ♦ ' MAUTZ  ,  JOE ' »MSGLEVEL=i 

//  EXEC  FORTGCIG , PARM.FORT= 1  MAP  1 
/ /FORT .SYS I N  DO  * 

COMPLEX  A3, A4,Z ( 1600) ,GS(40) ,G( 5292 ) ,U 
DIMENSION  RH(43) ,ZH( 43) ,DH(43) ,TJt20) 

DIMENSION  SV ( 42 ) , CV ( 42 ) , Z S ( 42 ) , R ( 42 ) , ANG ( 40 ) , AC ( 40 ) ,  CSM 120) 
DIMENSION  TP(80) ,T(80) ,TR(80) ,JK(4) 

REWIND  6 
U=(0,,1.) 

50  READ(l«51«EN0»b?)  NN,NP,NPHI ,  BK 

51  F0RMAT(3I3,E14.7) 

REAIH  1,53) IRH( I ) , 1  =  1, NP) 

RE AD ( 1,53) ( Z  H ( I ) » I =1 *NP ) 

53  FORMAT ( 10F8.4) 

76  WRITEt3,54)  NN ,NP .NPHI ,BK 

54  FORMAT ( IX//'  NN=',I3,'  NP=',I3,'  NPHJ=*,I3,«  BK=*  »£14,7 ) 

55  FORMATUX/'  RH') 

WRITE { 3, 55  ) 

WRITE ( 3 ,46 ) ( RH { I >,I=1,NP) 

46  r0RMAT(lX,10F8.4) 

WRITE ( 3  *56 ) 

56  FORMATUX/'  ZH'  ) 

WR I TE I  3 ,46 ) { ZH ( I ) , 1  =  1, NP) 

IF ( (RH( l }-RH(NP )).NE.O».OR«(ZH(1)-ZH(NP)).NE.O.)  GO  TO  58 
RH(NP+1)=P.H(2) 

ZH( NP+ 1 ) =  ZH( 2 ) 

RH (NP+2 ) =RH  (  3 ) 

ZH ( NP+2 ) =ZH { 3 ) 

NP=NP+2 

58  DO  57  1=2, NP 
12=1-1 

RR1=RH(I )-RH( 12) 

RR2=ZH(I  )-ZHU2) 

DH{  I2)=SQRT(RR1*RR1+RP.2*RR2) 

Z S t I2)=.5*(ZH( I )+ZH( 12)  > 

R( I 2)=.5*(RH( I )+RH( 12)  ) 

SV(I2j=RRl/DHU?) 

CV( 1 2)=RR2/0H(  125 

57  CONTINUE 
KG=NP— 1 
N=KG/2 
NM=N-1 
N«2=NM*2 
NM4=NM*4 
NZ=NM2#NM2 
NG=KG*KG 
M5=NN+2 
M6=NN+4 
FM=NN 
FM2=NN*NN 
PI=3. 141593 
ETA=376.707 
DP=P I /NPH I 
CA  =  BK*BK*ET  A 
CO=ETA 
SS=0. 

DO  117  1=1, NM 
Il=2#( I-1I+1 
12=11+1 

SS=$S+DH ( 1 1 ) +DH (12) 


11 


where 


COMPLEX  Z (NZ) 

DIMENSION  RH(NP) ,  ZH(NP),  U(NZ),  R(NZ),  T2INZ), 
A22 (NZ),  B(NZ),  X(NZ),  A(NZ),  Y(NZ),  T3(NZ), 
FKNZl,  EU  (NM2 ) ,  RU(NM2),  AMD  'NM2 ) ,  LB(NM2 ) , 

MB (MM2) 

NM2  =  NP  -  3 
NZ  =  NM2  *  NM? 


The  above  allocations  ax-e  based  upon  the  value  of  NP  after  execution  of 
statement  145 
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DO  13  K  =  1 ,NPHI 
K2=K+M4 

G!M2)=G(M2)+GS(K)*CSM(K2> 

13  CONTINUE 
68  CONTINUF 
17  CONTINUE 
16  CONTINUE 

DO  7 A  J=1 »NM 
J2  =  2  * ( J-l )  +  l 
J3= J2+ 1 
J4= J3+1 
J5= J4+ 1 
J6=4*< J-l H 1 
J7= J6+ 1 
J  0  =  J  7  + 1 
J9= J8+ 1 

DEL 1=DH ( J2 ) +DH ( J3 ) 

DEL2  =  DH(  J41+DHU5) 

TP(  J6)=DH(  J2I/DEI.1 

TP ( J7)=DH( J3I/DEL1 

TP ( J8)=-DH< J4I/DEL2 

TP( J9)=-DH( J5I/DEL2 

T( J6I=DH( J2)*DH( J21/2./DEL1 

T(J7)=DH(J3)*(DH(J2)+DH<J3)/2.»/DEU 

T  (  J8  •;  =DH  (  J4 )  *  ( DH  (  J5 )  +DH  <  J4  ) /?. .  I  /DEL2 

T( J9)=DH( J5)*DH{ J5I/2./DEL2 

74  CONTINUE 

DO  75  J=i,NM4 
TR( J)=T( J) 

75  CONTINUE 

115  I F ( (ZH( l)-ZH(NP-2) > .EQ.O. .AND. (RH( l)-RH(NP-2) ) .FQ.O. )  GO  TO  78 
IF(RH(1I)  77,23,77 

77  DEL1=DH ( 1 ) +DH ( 2 ) 

TR(  1)=DHI1  )*{  l.  +  (DH(2)+DH(  l)/2.  ) /DEL  15 
TR(2)=DH(2)*( l.+DH(2)/2./DELl ) 

23  IFIRH(NP) )  79,78,79 
79  Jl=  (  N-2  )  #4+3 
J2= Jl+1 

DEL2=DH (NP-2 ) +DH ( KG ) 

TR { J1 } =DH( NP-2 ) * ( 1 . +DH ( NP-2 ) /2./OEL2 ) 

116  TR I J2 ) =DH { KG ) # ( 1 .+ ( DHINP-2 ) +DH ( KG) /2 . I/D6L2) 

78  DO  30  J=1,NM 
JL=( J-1»#NM2 
J3=  ( J- 1 ) *4 

J  1=2=5=  ( j-l ) 

DO  31  1=1, NM 

L1=JL+I 

L2=Li+NM 

L3=NM*NM2+L1 

L4=L3+NM 

Z(L1)=0. 

Z (L2)=0. 

Z (L3)=0. 

Z (L4)=0. 

1 1=2=5=  ( I -1 ) 

1 3= ( I -1 ) *4 
DO  70  2 J= 1 ,4 
J2=J1+JJ 
J7=J3+J J 
DO  71  11=1,4 
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12=11+11 
J  7  “  I  3  *■  I  I 

j4=(  j?~ikkg+i? 

J5r J4+NG 
J  J  5+AJ6 
SS  =  SV ( I  2 ) *  SV ( J  ?  ) 

rc=cvi  m*cvu?) 

A3=.b*lG( JM+G(  J4)  ) 

A4=.5#(G<  JA)-G<  J4 )  ) 

7 (L 1 ) -Z ( LI )+*(CA*T (17) *T { J7  )  * ( $S* A3+CC*G { J5 ) )-CO*TP( I  7 ) «TP { J7 ) #6 ( J5 
l)KU 

7  U2)  =  Z { L2)+CA*SV( J2)*TR ( r 7 ) *T ( J7 ) *A4-FM*C0*G ( J5 ) *TR ( 1 7 ) *TP ( J7 ) /R { 
II?) 

Z  (  L  3  )  -  Z  (  L3)-CA*SV(  I  2 )  *T  U  7  )  *TR  (J7  )  *A4+FM*C0*G  (  J5  )  *-TP  {  1  7  )  *TR  ( J7  )  /R  { 
1J2) 

7  (  L  A  )  -Z  (  L4  )  +  (  C  A*A3-  FM2*C0/R  (I  2 ) /R ( J2 ) *G( J5 ) ) *TR  < 1 7 )  *TR  ( J7 ) *U 
71  CONTINUE 
70  CONTINUE 
31  CONTINUE 
30  CONTINUE 

WRITEIM  12  (  I  )  1 1  =  1  tNZ  ) 

88  FORMAT ( 1 X  » 10G11  .4  ) 

JK ( 1 ) =  L 
JK( 2)=N 
JK(3)=NM?*NM+1 
JK ( 4 ) = JK { 3 ) +NM 
DO  93  J  =  l,4 
K 1  =  JK ( J  ) 

WRITE (3, 24)  J 
24  FORMAT I  IX/*  ZM1) 

DO  92  I=1,NM 
K2=Kl+NM-l 

96  WRIT6(3  » 88  H  Z ( K) tK=Kl »K2) 

KUKI+NM2 

92  CONTINUE 

93  CONTINUE 
GO  TO  50 

52  STOP 
END 

/# 

//GO .FT 06F001  DD  0SNAME=EE0034 .REV1 f D I SP=OLD ,UN I T  =  23 14 1  X 

//  VOLUME  =  SER  =  SU0004  ?  DCB= ( RECFM^V , ELKS  I ZE= lHOOt  LRECL  =  1796 1 

//GO.SYSIN  DD  * 

1  21  20  0*3141 593E+00 
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III,  EIGENCURRENTS 

Program  112 ,  which  calculates  the  eigencur rents,  accepts  punched  card 
data  according  to 

READ (1,7)  NN,  NP 
7  FORMAT (213) 

READ  (1,39) (RK(I),  1  =  1,  NP) 

READ (1,39) (ZH(I),  1=1,  NP) 

39  FORMAT (10F8. 4) 

The  variables  NN,  NP,  RH,  and  ZR  have  the  same  meaning  as  in  program  111. 
The  generalized  impedance  matrix  Z  is  read  from  the  first  record  of  direct 
access  data  set  6  according  to 


REWIND  6 

READ (6) (Z (I),  1=1,  NZ) 

2  2 

where  NZ  is  either  (NP-1)  or  (NP-3)  depending  on  whether  or  not  C  closes 
upon  itself.  Program  ill  writes  the  eigencurrents  on  record  2  of  direct 
access  data  set  6. 

The  R  and  X  in  DO  loop  11  are  the  matrices  [R]  and  [X]  appearing  in 
(2-25)  of  [1]  computed  from  [Zn]  of  (2-47)  of  [1],  The  matrix  [X]  is  stored 
columnwise,  but  the  matrix  [R]  is  stored  according  to  the  symmetric  mode  used 
in  the  IBM  System/360  Scientific  Subroutine  Package  [4].  Statement  130  invokes 
the  eigenvalue  and  eigenvector  subroutine  EIGEN  in  the  Scientific  Subroutine 
Package.  The  subroutine  EIGEN  puts  the  eigenvalues  y  of  (2-26)  of  [1]  ordered 
y^  >_  >.  •  • »  in  the  diagonal  positions  of  R.  The  eigenvactors  of  the  matrix 

[R]  will  appear  in  U,  It  has  been  observed  that  the  eigenvectors  obtained 
from  EIGEN  are  normalized  so  that  U  is  orthogonal  in  accordance  with  (2-26) 
of  [1],  DO  loop  104  stores  the  y^  in  EU.  DO  loop  75  puts  the  matrix  [XU] 
appearing  in  (2-30)  of  [1]  in  T2.  DO  loop  78  puts  the  matrix  [A]  =  [U  XU]  of 
(2-30)  of  [1]  in  A,  Upon  exit  from  DO  loop  70,  JM  is  the  dimension  of  the  sub¬ 
matrix  [yp ]  in  (2-28)  of  [1],  DO  loop  73  puts  the  submatrix  [ A^2 3  appearing  in 
(2-30)  of  [1]  in  A^2»  Since  MINV  is  the  matrix  inversion  subroutine  from  the 
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Scientific  Subroutine  Package  14],  statement  128  inverts  the  matrix  [A  ]. 

-1  -1  ~ 

DO  loop  81  puts  t A22  ^  3  *n  ^3.  DO  loop  84  puts  [A^2][A22]  [A^J  ;‘-n 

B,  using  the  symmetric  mode  of  storage.  Statement  129  finds  the  eigenvalues 

and  eigenvectors  of  the  matrix  [B]  appearing  in  (2-36)  of  [1],  DO  loop  107 

-1/2 

puts  the  eigenvalues  A  of  (2-36)  of  [1]  in  AMD.  DO  loop  91  puts  y]  of 

(2-37)  of  [1]  in  T2.  Upon  exit  from  DO  loop  93,  the  matrix 

1«)  1  [w~i/2y] 

[a22  a12] 

will  be  in  T2. 

t* 

The  index  J  of  DO  loop  96  indicates  the  J  eigencurrent .  DO  loop  98 
puts  the  [I]  of  (2-37)  of  [1]  in  FI.  Because  f^(t)  of  (2-42)  of  [1]  is  de¬ 
fined  by  (30)  of  [3],  statements  143  and  144  have  to  divide  the  elements  of 
[I]  by  p  in  order  to  obtain  the  sinusoidal  components  of  the  eigencurrents. 
DO  loop  137  sets  the  largest  of  these  sinusoidal  components  equal  to  unity. 
At  the  time  FI  is  printed,  the  eigencurrent  (current  per  unit  length)  at 

t*  li 

the(2*J+l)  data  point  is  given  by 

t  FI  (J)  +  u,  FI (J+NM) 
t 

u.  FI(J)  cos  n$  +  u,  FI  (J+NM)  sin  n<{> 

c  9 

When  NN  =  n  ^  0,  the  alternate  eigencurrent 

u^  FI(J)  sin  n<£  -  u^  FI(J  +  NM)  cos  n<{> 
is  also  possible. 

If  NP  >  41,  some  dimension  statements  must  be  changed.  Minimum  allo¬ 
cations  are  given  by 


NN  =  0 
NN  =  n  j  0 


where 


COMPLEX  Z(NZ) 

DIMENSION  RH(NP),  ZH(NP), 
A22 (NZ),  B(NZ),  X(NZ), 
FI(NZ),  EU (NM2 ) ,  RU (NM2 ) , 
MBINM2) 

NM2  =  NP  -  3 
NZ  =  NM2  *  KM? 


U (NZ ) ,  R(NZ) ,  T2(NZ), 
A  (NZ ) ,  Y (NZ) ,  T3(NZ), 
AMD (NM2 ) ,  LB(NM2), 


The  above  allocations  are  based  upon  the  value  of  NP  after  execution  of 
statement  145 
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Listing  of  Proliant  • >2 

//  { 0034 <  £  E  *  2  *  2 ) » 'MAUTZ , JOE* ,MSGLEVEL=1 

//  LXbC  SSPCIG , P ARM , FOR  T  = 1  MAP ' 

/ /FOR  1  .SYS  IN  DU  * 

COMPLEX  ?  (  1444 )  ,Ut  ,U2,U3 

01  MENS  I  ON  PH (41) ,?H( 41 ) ,U( 1444) ,R ( 1444) ,T2( 1444) ,A22( 1444) ,B( 1444) 
0 1  MMS  I  ON  X(  1444)  » A  (  1444)  ,Y(  1444  )  ,  T3  ( 1444  )  ,  F  I  (  1444)  ,  EU  (  38  ) » RU  (  34  ) 
DIMENSION  AMD  (  3fl  '  « I.R  (  3ft  )  ,MB  (  3fi  ) 

EOH1 VALENCE  (R(l)»T2(l)*A22(l),8(l))t(X(l),A{l)»Y(in 
EOOI VALE Vf>  (TB( l ) ,P1 ( 1» ) , (EU( 1 ) ,AMD(1) ) 

PE  AIM  1,7)  MM ,  MP 
7  FORMAT  (ZH) 

WR  !  Tb  (  ? ,  3  )  NN,NP 
3  FORMAT ( *  1  NN  N P  * / 1 X , 2  I  3 ) 

READ!  1,39)  (RHU)  ,  1  =  1, NP) 

RFAD1 1,30) (ZH( I ) , 1=1 ,NP) 

39  FORMAT! iOPR. 4) 

HR  I T  E ( 3  » 4  0 ) ( RH ( I ) ,1  =  1 ,MP) 

40  FORMAT ( *  ORH  */(lX»10F8.4) ) 

WRITE (3,41) ( ZH ( I ) , I = 1 , NP ) 

41  FORMAT! *  07  H • / ( 1X,10F8.4) > 

145  IF! (RH! l)-RH(NP) ) .EG .0 . . AND . ! ZH ( 1 ) -ZH ( NP ) ) .EQ . 0 . 1  NP=NP+2 
144  NM2=NP-3 
NM=NM2/2 
NZ=NM2#NM2 
REWIND  5 

READ! 6) (Z( I ) , I=1,M? ) 

U3= ( 0 . , 1 . ) 

S 1  = .  2  5 

IF(MN.EO.O)  Sl= . 5 

U2=S1*U3 

J5=0 

DO  11  J  = 1 , NM2 
J?=( J-l )*MM2 
DO  12  I = 1 , J 
J5=J5+1 
J3= J2+ I 

J4= ( 1-1)*NM2+J 

IF(J.GT.NM.AND.I.LE.NM)  GO  TO  28 
U1=S1*(Z!J3)+Z(J4)) 

GO  TO  29 

28  Ul=U2*  (  7  ( .14  )  -Z  (J3)  ) 

29  R ( J5 1 =01 

X ! J3 ) = A I  MAG ( U 1 ) 

X  (  J4 )  =  X  ( ,J3  ) 

1?  CONTINUE 
11  CONTINUE 

130  CALL  E I  GEN! R ,U  »NM2  *0 ) 

J1=0 

00  104  J= ] , NM2 

J1=J1+J 

FU( J ) "R ( J1  / 

RU!  J)  =  l./SORT(  ABSIEIJ!  J)  )  ) 

104  CONTINUE 

WR I TE ! 3 , 14 1 ) ( EU ( J ) , J= 1 , NM2 ) 

141  FORMAT! 'OEIGENVALUES  OF  THE  MATRIX  R ' / ( IX  ,  7E 1 1 .4 ) ) 

DO  75  J'-ltNM? 

J 1 = ! J~l)#NM2 
DO  74  I = 1 ,NM? 

J?  =  J1-U 
T2( J2)=0. 


13 


»J3=  (  1-1 )*NM2 
OH  7?  K-lfNM? 

K 1  =K  + J3 
K2=K+ J 1 

T2(  J2)=T2(  J2)+X(K1)*U(K2) 

77  CONTINUE 
76  CONTINUE 
75  CONTINUE 

00  78  J“1 » NM? 

Ji=( J~l)*NM2 
OH  79  1  =  1  f J 
J2«J1+I 
A( J2)=0. 

J  3  = ( l-l )*NM2 
00  80  K=1,NM2 
K 1  =  K  + J3 
K?=K+J1 

A(  J2)=AU?)+U(K1  )«T2(K2) 

80  CONTINUE 
J4= J3+J 

A(  J4)=A( J?) 

79  CONTINUE 

78  CONTINUE 
X2=EU( l)«l*E-03 
DO  70  jsi9NM2 
JM= J—l 

IE(EIJ{J)  .LT.X2)  GO  TO  ?2 
70  CONTINUE 

72  JN=NN2-JM 
JM I = JM+ 1 
J1*0 

00  73  J=JM1,NN? 

( J-1 )*NM? 

00  74 

J1=J1+1 

J3= J2+ I 

A22 ( J 1 ) = A ( J3  ) 

74  CONTINUE 

73  CONTINUE 

I?8  CALL  MINV(A22t JNtl)«LB9MR) 

J1  =  0 

no  81  J=1,JM 
J3=( J-1 )*NM2+JM 
00  82  1=1 tJN 
J?  =  t 1-1 )*JN 
J1=J1+1 
T  3  (  J  1 )  =  0  * 
no  83  K=  1 1  JN 
K  1=J2+K 
K2  =  J3*rK 

T  3 ( J 1 ) =  T 3 ( J 1 )+A?2(Kl)#A(K2) 
83  CONTINUE 
8?  CONTINUE 

81  CONTINUE 
J?  =  0 

no  84  J=1,JM 
J3= ( J-1 ) #NM2 
J5~ ( J-1 ) *JN 
DO  85  1=1, J 
J2= J2+ 1 


J6=J3+I 
H  ( ,12  )  =A(  J6  ) 

J6=(  1-1  )  #Mi '2+.IM 
lltl  A  A  K=1,JM 
K  1  =  K  +  J  6 
K2=K+J6 

3(  ,!?)=(<(  J?)-A(K1  ) *1 3 (K2  < 

36  CONTINUE 

B(J?)=B(  .12)  MMH  J  }>.<RU(  I  ) 

86  CMNTIMUF 
86  CONTINUE 

1?9  CALL  FIGFNIB.Y. JM,0) 

J1  =  0 

(Ml  107  J“1  *  JH 

Jl  =  Jl+,l 

AMI )  (  J  )  =  8  (  J  1  ) 

107  CONTINUE 

WRITE (3* 58)  (  AMO  (  J  )  ,  J=1,JM) 

58  FORMAT! 'OF  IGr.NVALUES  OF  THE  MATRIX  B  '  /  (  IX  ,  5E  16 . 7 )  ) 
00  91  J=  1  ,,IM 

ji=i j-i )*jm 

.16=  (  J-l  )*NM2 
00  9?  1=1  ,.IM 
.13  =  I  +  J6 
J  2  =  I  +  J 1 

T  2 ( J3 ) =Y ( J?  )  *RU (  1  ) 

9?  CONTINUE 
91  CONTINUE 
S1=0. 

on  93  J=l, JH 
J 1 = ( J-1)*NM2 
on  96  I  =  1 . JN 
J2=J1+I +JM 
121.121  =  0. 

00  9  5  K  =  1 .  JM 
K  1  =  (  K- 1  )#.IN+I 
K2=K+J1 

12  1 J2)=T?( J21-T31K1 )*T2(K2) 

95  CONTINUE 

96  CONTINUE 
93  CONTINUE 

011  9 A  J=1.JM 
S1=0. 

Jl=( J-l)*MM? 

DU  97  1=1 »NM 
,I2  =  J1  +  I 
J3= J2+NM 
,16  =  2*1  +  1 
FI  ( J?)=0. 

FI  1 J3)=0. 
on  98  K=1,NM2 
K2  =  K+.I 1 
K 1  =  ( K-l ) "NR  2  + 1 
K3=K1+NM 

FI ( J2)=F1 1 J21+UCK1 )#T2(K2) 

F! ( J3)=FI ( J3)+U(K3)*T2(K2) 

98  CONTINUE 

163  FI  (  J2)=F1  {  J?)/P.H(  J6) 

166  El (J3)=FI ( J3)/RH( J6) 

S2  =  APS  (  F  1 1 .12 )  ) 
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S3=Aft$ ( F I ( J3)  ) 

I F ( S2-S1 )  136,133,133 
133  Si=S2 
J5=J2 

136  1 F ( S3-S 1 )  97 , 135 • l 3b 
13b  S l =S3 

J  5 = J  3 

97  CUNT  I NUF 

Sl=l •/FI ( Jb) 

J2=J1+1 
J3sJl+NM 
=  J3+ 1 
J5= J 1+MM2 

WR I TF ( 3 t 138)  AMD(J) 

138  FORMAT  v  1 OL AMBDA  =  Sbll.4) 
nn  137  I  =  J  2  ,  J  5 
FI (I )=FI ( I  )*S1 

137  CI1NTIMUE 

WR I TE ( 3 , 60 ) ( F I ( I  ),I=J2,J3) 

60  FORMAT ( 1  JT *, 1 0F8 . 4 / ( 3X , 10F8 ,4 ) ) 

WR  I  T  E  (  3 , 6 1  ) 

61  FORMAT (*  JO') 

WRITE (3,62) (FI (  I  )  ,  I=J4,J5) 

62  FORMAT ( * +  / 1 , 10F8 * A/ { 3X , I 0F8 .4 ) ) 

96  CONTINUE 

WRITE (6) (FI ( I ) • 1=1, J5) 

STOP 

END 

/# 

/ /GO » FT O6F00 1  OD  D$NAME=EE0034 . REV1 , DI SP“OLO, UN  I T=23 14 ,  X 

//  VOLUME2  SER  =  SU0004~,  DGB=  (RECFM=V , BLKSI 7E=1B00 ,LRECL~i796) 

//GO. SYS  IN  00  * 

1  ?  1 


n.oooo 

0.5000 

1.0000 

1.5000 

2.0000 

2.5000 

3.0000 

3.5000 

4.0000 

4.5000 

5.0000 

10.0000 

5.5000 

6.0000 

6.5000 

7.0000 

7.5000 

8.0000 

ft. 5000 

9.0000 

9.5000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

/* 

// 


6.or;)j  6.^0.*^  7  *0j0)  7,s^o 

iC«u*K  *> 
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Output  oi  Program  f*2 
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IV,  EIGENCURRENT  PLOTS 

Program  //3  which  plots  the  eigencurrents  accepts  punched  card  data 
according  to 

READ (1,90)  NM2,  JM 
90  FORMAT (2 013) 

In  the  previous  program  an  NM2  by  NM2  generalized  impedance  matrix  has  led  to 
JM  eigencurrents.  More  precisely,  NM2  is  defined  by  statement  146  of  the 
previous  program  and  JM  is  the  number  of  numbers  printed  under  the  heading 
"Eigenvalues  of  the  matrix  B"  in  the  output  of  the  previous  program.  The 
eigencurrents  are  read  from  the  second  record  of  direct  access  data  s<st  6 
according  to 

REWIND  6 
READ  (6) 

NZ1  =  NM2  *  JM 

READ  (6) (FI  (I),  1  =  1,  NZ1) 

DO  loop  93  prepares  the  vertical  coordinate  for  plotting  and  DO  loop  95 
prepares  the  horizontal  coordinate.  The  origin  is  at  (1,,5,),  The  hori¬ 
zontal  axis  corresponds  to  the  contour  length  variable.  If  the  NP  data 
points  defining  the  contour  C  are  not  equally  spaced,  this  correspondence 
may  become  nonlinear  because  DO  loop  95  always  supplies  equally  spaced 
horizontal  coordinates.  The  index  J  of  DO  loop  94  indicates  the  J 
eigencurrent .  Statements  18  and  19  plot  the  sinusoidal  components  of 
and  respectively. 

Minimum  allocations  for  FI  and  X  are  given  by 


DIMENSION  FI (NM2  *  NM2),  X(NM2/2). 
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Li&Liuu  of  l'nif.rui  M 

//  ( 003a  «  FF * MAUTZ  ,  JME  *  ,MSG1JEVEL  =  1 

//  EXEC  FORT  GCLG  »  PARim  ,  FUR  T  =  1  MAP  • 

//FOR  I  .  SY -s  I  m  0U  : 

DIMENSION  I-  1  (  1  AAA  )  ,XP  ( A  )  ,YP(A  )  ,AREA< A00)  ,X(  19) 
f.ALl  PLOTS!  AOrAfAf*OJ 
READ!  1  ,*>(!) 

40  FORMAT  (  ?0I  3) 

WRITE  <3  ,01) 

4!  FORMAT  I  U>NM?  JO,' /IX,  21  3) 

0/1=00?  UM 
PFWlNO  6 
REAM  M 

RFADI61  (  F  I  ( I  ) ,  !=l,N71 ) 

FM=MM?/2 
XP( 11=1* 

XPI21=6. 

YPI  l)=t>. 

YP(2)=S. 

XP( 31  =  1  . 

X  P  ( A  1  -  1  . 

Y  P  (  3  )  =  3  . 

YP(a>=7. 

DO  43  !  =  1.M71 
HI  I  )=5.+El <  I  ) 

43  CUNT  I NUF 

Sl  =  t>./(NM*1) 

UO  95  J=i,NM 
X | J ) = 1 . + J*S 1 

95  CONTINUE 

I/O  9 A  ,1=)  ,  JM 
Jl=( J-l)*NM2+l 
j?=ji+riM 

CALI.  LINE (XP(l) ,YP(1) *2 *1,0,0) 

1)0  96  K  =  ?,6 
S1=b-k 

CALL  SYMBOL  (  51 ,5  .  .  . 1 A . ) 3 ,0 . ,-l > 

96  CONTINUE 

CALI  LlfcfcUP<3),Y«M3),2,l,0,0) 

00  97  K=  1 , 5 
S1=B-K 

CALL  SYMBOL II., SI , . 1A , 13,90. ,-T , 

97  COM  I Nl IF 

CALL  NUMBER ( .76,5.93, .1A,1.,0.,-1 i 
CALL  NLHBERI  .76,4,93,  .1A,0.,0.,— 1  ) 

CALL  NUMBER ( .6 A, 3. 93, .14,-1. ,0.,-l ) 
lft  CALL  LINPIXII) ,FI I Jl> ,NM,1,4,1) 

19  CALL  LIMF.IXII)  ,FI( J?)  ,NM,  1,0,1) 

CALL  PLOT  1 7., O., -3) 

4<,  C  OMT  I  O'  IF 

C&LL  PLOT  I  5,, 0.. -3) 

STOP 

ENO 

/v 

//GO.ETO6FOO1  00  DSNAM6  =  EE003A .RFV1 ,D 1 SP=OLD , UN  I T=23 1A ,  X 

//  VOLUME  =$ER  =  $U0004, UCflMRECPM=V,BLKS  I ZE=lflOO,LRECL  =  l 796) 

//UO.SYSIN  on  -j 
U  3 
/* 

// 
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V .  GAIN  EIGENPATTERNS 

Program  //4  which  calculates  and  plots  the  gain  patterns  of  the  eigencur- 
rents  accepts  punched  card  data  according  to 

READ (1,10)  NN,  NP,  NT,  NS,  JM,  BK,  SCL 

10  FORMAT  (513,  2E14.7) 

READ  (1,11) (AMD (I),  1  =  1,  JK) 

11  FORMAT  (5E14.7) 

READ (1,15) (RH(I),  I  =  1,  NP) 

READ(1,15) (ZH(1),  1=3,  NP) 

15  FORMAT  (10F8.4) 

Here,  NN,  NP,  BK,  RH,  and  ZH  are  the  same  as  in  program  //I  and  JM  is  the 
same  as  in  program  it 3.  The  electric  field  and  gain  will  be  computed  at 
polar  angles  0.  =  i=l , 2 , . . .  NT  but  will  be  printed  only  at  i=l, 

NS  +  1,  2*NS  +1,...  .  One  inch  will  correspond  to  a  gain  of  1./5CL  on 
the  plot.  The  variable  AMD  appears  under  the  heading  "Eigenvalues  of  the 
Matrix  B"  in  the  output  of  program  //2.  The  JM  eigencurrents  are  read  from 
record  2  of  direct  access  data  set  6  according  to 

REWIND  6 
READ (6) 

READ (6) (FI (J),  J  =  1,  NZ1) 

where  NZ1  -•  NM2  *  JM.  The  variable  NM2  is  the  same  as  that  appearing  in  the 
first  three  programs. 

In  DO  loop  40,  DH(I),  RS(I),  ZS(I),  SV(I),  ‘and  CV(I)  are  respectively 

t*  Vi  j-u 

the  distance  between  the  I  n  and  (1+1)  data  point,  and  p,z,  sin  v,  and 

cos  v  midway  between  the  I1"*1  and  (I+l)^'  data  points.  Here,  v  is  the  angle 
between  u  and  the  z  axis.  In  DO  loop  41,  T(4*(J-1)  +  I),  1=1, 2, 3, 4,  is  the 

Tj(t) 

amplitude  of  -  (the  J  triangle  function  divided  by  the  cylindrical 

coordinate  p)  multiplied  by  DH(2* (J-l)+I)  and  evaluated  midway  between  the 
(2*(J“1)+I)fck  and  (2* (J-lJ+I+l)*"^  data  points.  Now  T  represents  the  func¬ 
tions  f^(t)  associated  with  J^  and  TR  those  functions  £^(t)  associated  with 
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although  only  f^(t)  appears  in  (2-42)  of  [1],  no  distinction  being  made 
there  between  f*T(t)  and  JL't(t).  If  the  surface  S  has  no  edges,  f^(t)  =  ff(t)r 
but  if  S  has  an  edge  at  either  end  of  C,  program  //4  modifies  the  f?(t)  nearest 

i  2  i 

this  edge  so  that  ~  at  the  edge  in  an  attempt  to  account  for  the  singu¬ 

larity  of  J  at  the  edge. 


The  subroutine  PLANE  is  concerned  with  the  integrals 

jk  *  r 

^  3  "i  *  \  G  3  m  ds 


appearing  in  (2-51)  of  [1],  The  R  integrals  can  be  written  as 


<«- 


The  row  submatrices  on  the  right  hand  side  of  (3)  are  defined  according  to 


which  W.  and  which  u  are  used 
i  m 


J  Matrix  element 

rt  t6 


~y 

u 

m 

-y 

w. 

1 

ue 

—y 

u(.fi(t)  cos  n<j> 

->- 

ue 

u^f^(t)  sin 

-y 

•> 

U<1> 

u^f i  (t )  cos  n(j» 

-> 

-> 

U, 

<P 

u^fj.Ct)  sin  n<j> 

If  the  coefficients  appearing  in  (2-51)  of  [1]  are  partitioned  into  column 
vectors  IC  and  1^  corresponding  to  the  ut  and  u^  directed  W^,  the  radiation 
field  will  be  given  by 
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* 


-ja)ue 


-j  kr 


4nr 


c 

c 

c 

Furthermore,  it  will  be  shown  that 
E, 


e 


lE4 


-jojue 


-j  kr 


47Tr 


cos  n4  0 

0  sin  n4 


[Ifc] 

[1^] 

(4) 


ite 


,  ji6* 

L  J  n 


-jR*6 
J  n 


R 


4*6 

n  J 


I1 

M 

rr 

- J 

1 

M 

-e- 

4  -.1 

1 

(5) 


where  the  careted  submatrices  R^  are  given  by  equations  (77)  and  (81)  of 

reference  [3].  If  the  current  has  the  polarization  (2-43)  of  [1]  instead 

of  (2-42)  of  [1],  the  matrix  [C0^  of  (5)  must  be  replaced  by 

[Sin  The  submatrices  use  sinusoidal  VL  and  depend  upon  the 

measurement  azimuthal  angle  4  while  the  R  use  exponential  functions  and  are 

m  n 

evaluated  at  4  =  0.  Consider  the  contribution  to  En  from  IT,  According  to 

m  u  i 

(4), 

-j  kr  r  ^ 

E0  =  -  4  J  p  dt  j  d4-  ^ (9)^ (t)  sin  n<{>  •  (4>-<J>m)  (6) 

0 


-y  -y 

where  -j  k  •  r 

=  e  m  (7) 

Since  the  integrand  of  (6)  is  periodic  in  4  with  period  2^r,  4^  can  be 
added  to  4  without  changing  the  value  of  the  integral. 


-j  kr  2?r 

\  1 4  \«  d*  «t (♦*+„)  •  ;,<♦»>  ♦<*) 


sin  n(4+4m)  (8) 


But 


"*<♦+♦■>  •  VV  =  V*}  *  “e(0) 


(9) 


so  that 
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E  =  zte — 

af»  4nr 


p  f i (t )  dt  J  dd>  u  (♦)  sin  n (4>+^m)  *  u0(O)  lK$) 


Now  Up  (0)  is  the  measurement  plane  wave  coming  from  the  direction  <j>  =  0. 

Alternatively,  (10)  can  be  obtained  directly  from  (6)  by  reasoning  that 

■> 

(6)  should  depend  only  on  the  phase  difference  between  the  current  J  and 

the  direction  6  from  which  the  measure  ont  plane  wave  comes*  For  instance, 
m 

one  should  be  able  to  turn  both  J  and  4>  back  by  $  to  obtain  (10),  Equa- 

m  m 

tion  (10)  leads  to 


.  -jkr  .  R+6  -  R*6  R+6  + 

■JSjsJ—  cos  n*  +  S—r 

4Ttr  i  4.j  m  i 


R^°  +  R*0 
n  -n 


sin  n (j> 


Using  the  fact  that  is  odd  in  n, 

E  =  —  if  (-j  R^9  cos  n<f>  )  1 

e  fax  i  J  n  m 

in  agreement  with  (5).  The  rest  of  (5)  can  be  similarly  verified. 

The  subroutine  PLANE  is  essentially  the  same  as  the  one  appearing  in 
Appendix  B  of  reference  [2].  For  the  Lth  measurement  polar  angle  6  •-  THR(L) 
and  the  JCh  function  fj(t)  of  (2-42),  PLANE  stores  R^  ,  -jR^  ,  jRn  ,  and  Rfl 
in  VVR(L1+J),  VVR(L1+J+NM) ,  VVR(L1+J+NM*2) ,  and  VV:<(L1+J+NM*3)  respectively 
where  LI  =  (L-1)*NM*4  and  NM  =  NM2/2. 

DO  loop  92  multiplies  the  eigencurrents  by  p  to  retrieve 


[I]  = 


appearing  in  (2-25)  of  [11.  The  Lth  elements  of  [TC]  and  £ T^]  are  the  coefficients 
of  the  expansion  functions  ’£^(0  cos  n<{s  and  u.f^t)  sin  n<}>  respectively  for 
an  eigencurrent  of  polarization  (2-42)  of  [1]. 
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til 

The  index  KK  of  DO  loop  81  indicates  the  KK  eigencurrent .  In  DO  loop 

83,  K=1  obtains  and  K=2  obtains  E^.  DO  loop  82  obtains  the  NT  polar  angles 

0.  Inner  DO  loop  84  performs  the  actual  matrix  multiplication  indicated  by 

2  2 

(5).  The  gains  and  are  proportional  to  |E^|  and  j  E^  j  .  The  average 
of  the  total  gain  over  the  area  of  the  radiation  sphere  is  unity. 


1 


(Ga(0)  +  G, (0))  sin  6d0  = 
0  © 


NN  =  0 
NN  >  1 


(14) 


DO  loop  85  stores  E^  and  G^  in  positions  1  to  NT  of  E  and  G  and  E^  and  G 

positions  NT+1  to  2*NT  of  E  and  G.  The  phase  of  E  is  normalized  to  -je 

2  2 

The  magnitude  of  E  is  normalized  so  that  | E^ |  =  G^  and  |e^|  =  G^.  Statements 

34  to  86  are  concerned  with  plotting  the  gains  of  the  eigencurreuts. 


Minimum  allocations  are  given  by 


COMPLEX  VR(NT*NM2*2) ,  E(NT*2) 

COMMON  RS (NP-1) , ZS (NP-1) ,  SV(NP-l) ,  CV(NP-l) ,  T(NM2*2) ,  TR(NM2*2) 
DIMENSION  AMD (JM) ,  RH(NP) ,  ZH(NP) ,  DH(NP-l) ,  TH(NT) ,  G(NT*2) ,  SN(NT*2) 
CS (NT*2) ,  FI (NM2*NM2) ,  GX(NT*2) .  GY (NT*2) 

DIMENSION  BJ(3*(NP-1) ) 

Here,  NP  is  the  value  of  NP  after  execution  of  statement  96  in  the  main 
program.  Also,  BJ  appears  in  PLANE. 


Listing  of  Progxam  04 

//  ( 0034  ,-£6,3,2*  .5) ,  '  MAUT7 ,  J0E'«  ,MSG1EVEL=1 

//  EXEC  FORTGCIG,PARM.FORT= 'MAP' 

//F0RT.SYS1N  DD  * 

SUBROUTINE  PI ANE ( VVR , THR , NT  ) 

COMPLEX  VVR( li  ,A5,A6,U 

COMMON  IJ,  RSI  40)  ,ZS(40)  ,SV(40)  ,CV(40)  , BK, NP ,NN , T ( 80 )  ,TR(80) 
DIMENSION  B J ( 1 ?6 ) , THR ( 1 ) , FK ( 20 ) 

KG=NP-1 

NM=KG/2-l 

M2=NN+2 

A5=2  .#3 . 14 lb93*U**(NN+l ) 

NV=NM*4 
r-K(i)  =  i. 
no  153  J=1,M2 
J1=J+1 

FK ( J 1 ) =FK(  J ) 
i 53  CONTINUE 

DO  156  L  =  1 » NT 
L 1= ( L-l ) #NV 
CS=COS ( THR ( L ) ) 

$N=S I N ( THR { L ) ) 

BCS=BK~CS 
DO  302  J=1 *KG 
X=RS( J)*BK*SN 
J1=J 
I  1  =NN 

1F1I1)  303,304,303 

304  11=11+1 

J  1  =  .J  1+KG 

303  DO  305  J J= 1 1 ,M2 
IFIX-l.E-5)  1,1,2 

1  IF(JJ-l)  3,3,4 

3  BJ ( J 1 ) =1 , 

GO  TO  306 

4  B.!(J1)=0* 

GO  TO  306 

2  RH=X/2. 

RH2=RH=f=RH 
RH3=RH** ( J J-l ) 

BJU1)=RH3/FK(  J  J ) 

SS=B J  t  J 1 ) 

R  SST  =  SS*=1 .6-7 
DO  155  K=1 ,20 
SS=-SS*RH2/K/(K+JJ~1) 

BJ( J1)=RJ( J1)*SS 
IF( ABS(SS)-SST)  306,306,155 
155  CONTINUE 
STOP  155 

306  J  1=J 1+KG 

305  CONTINUE 
302  CONTINUE 

I  FINN)  307,306,307 

308  DO  309  J=1,KG 
Jl=J+t*KG 

BJ ( J ) =-BJ ( J 1 ) 

309  CONTINUE 

307  DO  300  J=  1  »NM 
J 1  =  J+L 1 

J2= J 1+NM 
J3= J2+NM 
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J4= J3+NM 
VVR ( J 1 ) =0. 

VVR! J2)=0. 

VVR  (  J3  )  '-0 . 

VVR ( J4 ) =0 . 

DO  301  1=1,4 
1 1  =  2* ( J-l )  + 1 
14=4* ( J- 1 ) + I 
I  2=1 l+KG 
13=1 2  +  KG 

A6= ( COS ( Z  S ( I  1 >*RCS)+U*SIN(ZS( I 1)*BCS) ) *A5 
B  J 1  =  (  B  J  ( I  3 ) +RJ ( 1 1 ) ) *.  5 
b J2=(BJ( 1 3 )-BJ ( 1 1  )  )*.5 

VVR( J1)=VVR( J1 )+A6*(CS*SV( I 1 ) *BJ2+SN*CV ( 1 1 ) *BJ ( 1 2 ) *U ) *T ( 1 4 ) 

VVR!  J2)=VVR( J2)+A6*CS*BJ1*TR(  14) 

VVR  (  J3  )  =  VVR  ( .13  )  +A6*SV (  I  1  )*BJ1*T(  14) 

VVR  ( .14  )  =VVR  (  J4  )+A6*BJ2*TR  (  14  ) 

301  CONTINUE 
300  CONTINUE 
156  CONTINUE 
RETURN 
END 

COMPLEX  U  ,U  1  ,VR! 5548 ) ,  E ( 146) 

COMMON  U, RSI  40) ,ZS(40) ,SV(40) ,CV(40) , BK ,NP , NN , T ( 80 ) ,TR(80) 
DIMENSION  AMD (38 ) »RH(41 ) ,ZH<41) ,DH(40) ,TH(73) » XP ( 2  )  »YP(2)  ,G( 146) 
DIMENSION  SN ( 73 ) ,CS ( 73 ) , F I ( 1444 ), GX ( 146 !, GY ( 146 ), AREA ( 400 ) 

CALL  PLOT S ( AR EA , 400 ) 

READ! 1, 10)  NN,NP,NT,NS, JM,BK,SCL 

10  FORM AT i 5 1  3, 2E 14. 7) 

READ! 1, 11) (AMD! I ) , 1=1, JM) 

11  FORMAT  I 5E14.7 ) 

RF  ADI  1 , 15 ) ( RH ( I ) , I  =  1 , NP ) 

READ! 1,15) (ZH( I ) ,I=1,NP) 

15  FORMAT! 1CF8. 4) 

WR  ITE ( 3 ,33 ) 

33  FORMAT  I  *  1  NN  MP  NT  NS  JM » , 6X , ' BK • , 12X, ' SCL ' ) 

NR  ITE (3, 12)  NN  ,NP  »NT  »NS , JM  ,  BK  » SCL 

12  FORMAT! IX, 51 3 , 2E14.7 ) 

WRITE (3,5) I  AMD! I ) , I  =  1 , JM ) 

5  FORMAT! 'OAMD'/( 1X,5E14.7) ) 

WRITE(3,1A) (RH(I ) , I = 1 ,NP ) 

16  FORMAT! • ORH ' / (1 X , 10F8 . 4) ) 

WR  ITE (3 , 18 ) ( ZH ( I ) ,1  =  1 ,NP) 

18  FORMAT! » OZH 1 / I IX , 10F8 . 4) ) 

U=(0.,1.) 

PI=3. 141593 
KL  =  1 

I F ( (RH( l)-RH(NP) ) .NE.O..OR. ( ZH ( 1 )-ZH(NP) ) .NE.O. )  GO  TO  96 
KL  =0 

RH ( NP+1 ) =RH ( 2 ) 

ZH(NP+1)=ZH(?) 

RH ( NP+2 ) =RH ( 3 ) 

ZHINP+2)=7.H13) 

NP=NP+2 
96  NM2=NP-3 
NM=NM2/2 
NM4=NM*4 
NT2«NT*2 
NZ1=NM2#JM 
REWIND  6 


READ! 6) 

Rf-ADC  6)  ( F  I  (  J  >  .  J  =  )  ,  N  Z  1  ) 

DM  4-0  I  -?  »MP 
12=1-1 

RR1=RH(  D-RHII?) 

RR2=ZH( 1 )-ZH( I?) 

DH(  1  ?  )=SORT  (RR1*RR1+RR?*RR2  ) 

R  S ( 1 2 )  =  . 0* ( RH ( I ) +RH (12)) 

ZSU2)  =  .b*<ZH<  I  )+ZH(  I?)  ) 

SV( I2)=RRl/DH( 1?) 

CV( I?)=RR2/DH( 12) 

40  CONTINUE 

DC)  41  J  =  1,NM 
J2=2*( J-ll+1 
J3= J2+1 
J4=J3+1 
J5= J4+1 
J6=4* ( J- 1 ) + 1 
J7=J6+1 
J8= J7+  1 
J9= JR+1 

DEL 1=DH( J2 ) +DH ( J3 ) 

DEL2=DH( J4)+DH(Jb) 

T ( J6 ) =DH ( J2 )  *l)H  (J2J/2./DEL1 
T ( J  7 ) -DH ( J3 ) * ( OH (J2 ) +DH ( J3 ) /2 . ) /DEL  1 
T( J8)=DH{ J4)*(DH( J5)+DH{ J4)/2. )/DEL2 
T(  J9  )=DHUb)*DH(.l5)/2./DEL2 

41  CONTINUE 

DO  97  J=1,NM4 
TR ( J ) =T ( J ) 

97  CONTINUE 
IF(KL.EO.O)  GO  TO  98 
IF(RHll))  23,24,  23 

23  DELI =DH ( 1)+DH(2; 

TR ( 1 ) =DH ( 1 ) * (  l.+  (DH(2)+DH( l)/2.)/DELl) 
TR(2)=DH(2)*( l.+DH(2)/2./DELl) 

24  I F ( R H ( N P ) )  26,98,26 
26  Jl=(NM-l)*4+3 

J2= J 1+1 

DEL2=DH { NP-2 ) +DH( NP-1 ) 

TR( J1 )=DH(NP-2)*( 1 ,+DH(NP-2)/2./DEL2) 

TR( J2)=DH(MP-1)*( 1 .+(DH(MP-2)+0H(NP-l)/2. )/DEL2) 

98  DEL=P I / ( NT- 1 ) 

C2=SCL*4 . /DEL 
IF(NN.EO.O)  C 2=C?=» .5 
DO  43  J= 1 ,NT 

TH( J ) = ( J-l ) $DEL 
SN( J)=SCL*SIN(TH(  J) ) 

CSU)=SCL*COS(THU! ) 

43  CONTINUE 
XP ( 1 ) =2 ♦ 

,vp !  2 )  =  8 » 

YP '( 1 ')  =  5 » 

YP ( 2 ) =5 . 

CALL  PLANE ( VR ,TH ,NT ) 

C1=180./PI 
UO  31  J=1,NT 
TH{ J)=TH( J)«C1 
31  CONTINUE 

DO  92  J=1,JM 


Jl  =  ( J-1)*NM2 
00  93  1=1 »NM 
J2=J 1+1 
J3= J2+NM 
J4=2*I+1 

FI ( J2 ) =F I ( J  2 ) SRH ( J4 } 

FI  ( J3  )  =  FI  I  J3  )  ,!<RH  (  J4) 

93  CONTINUE 
9?  CONTINUE 

00  81  :.is=l,JM 
Kl= ( KK-1 ) *NM2 
S2  =  0. 

00  83  K  =  1 » 2 
K2= ( K— 1 ) *NT 
K3= ( K  —  1 ) *NM2 
00  82  J=1,NT 
J 1= ( J-1)*NM4+K3 
U1  =  0. 

00  84  1=1, NM2 
J3=I+J1 
J4=I +K1 

U1=U1+VR( J3)*FI ( J4) 

84  CONTINUE 
J5= J+K2 

E ( J5 ) =U1 
S1=CABSIU1) 

G(J5)=S1*S1 
S2=S2+G( J5)*SN( J) 

82  CONTINUE 

83  CONTINUE 
S3=S0ftT  (C2/S2) 

00  85  J=1,NT2 
E( J)=E( J)#S3 
S2=CABS(E( J) ) 

G I J ) =S2*S2 

85  CONTINUE 

34  CALL  LINEIXP, YP, 2, 1,0,0) 

00  77  J= 1  *  7 
Sl=9-J 

CALL  SYMROL (SI, 5. ,.14, 13,0. 

77  CONTINUE 

CALL  LINE(YPtXP, 2,1,0, 0) 

DO  78  J=1 ,7 
S 1=9-J 

CALL  SYMB0LI5. , SI ,  .14,13,90 

78  CONTINUE 

00  86  K=1 , 2 

K2= ( K-l ) *NT 

DO  87  J=1 »NT 

Jl= J+K2 

S1=G ( J1 )*SN5  J  5 

GX I J )=5.+Sl 

GY  I J )  =5  ,.+G  (  J 1 )  *CS  (  J ) 

J2=NT2-J+1 

GX( J2)=5.-S1 

GY ( J2 ) =GY ( J ) 

87  CONTINUE 

CALL  LINE(GX,GY,NT?,1,0,0) 

86  CONTINUE 
WRITE(3,27)  AMO ( KK ) 
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2  7  HiHMAI  <  ‘OFLtCTKIC  Hell'  AND  GAIN  OF  fc  I GENCURRFNT  FOR  LAMBDA  =',011 

1  .-) 

l-'K!  IK  3,2k) 

FOR  MAT  CO  0  Rf-AL(FO)  I  MAG  ( 80 )  GAINO  REAL  ( 00 )  I 

CAGIEO)  GAINO') 

WR I TF ( 3 , 29 ) 

29  FORMAT  (  '  +  "',10X,'-',11X,'-',11X,'-',11X,'/',11X,'/',11X,'/') 

no  JO  K=1,NT,MG 
8  2-N T +K 

W«  I  ri-  (  3  *  32 1  TH(K),F(K),G(K),F(K2),G(K2) 

32  FORMAT  (  1X,F6.1 
*.0  f.ONT  I NlfH 

CALL  PLOT ( 7 . ,0. «-3) 

«1  CONTI  NOE 

CALL  PLOT ( 6 . ,0, » -3 ) 

STOP 

FMIT 

/  1 

//GO.FTOhFOOl  DO  DSNAME  =  000034  ,RE VI ,01 SP=0L0 ,UN I T  =  23 14 ,  X 

//  VOlOM0=SER=ROOOO4,nCR=(RECFM=V,BLKSIZB=lfiOOf LRECL=1796) 

//GO. SYS  IN  DO  F 

1  21  73  4  3  0.3141 593B+00  0. 5000000F+00 

0.25596650+0 1-0.87856840-02-0.24445270+02 


0.0000 

0.5000 

1.0000 

1.5000 

2 .0000 

?  .5000 

3.0000 

3.5000 

4.0000 

4.5000 

6.0000 
10.  noon 

5.5000 

6.0000 

6.5000 

7.0000 

7.5000 

8 .0000 

8.5000 

9.0000 

9.5000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0„0000 

0.0000 

o.onoo 
O.Oi 00 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000, 

0.0000 

/* 

// 
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Output  of  Program  /M 

KH  NP  NT  NS  JM  H*  SCL 

L  21  73  4  3  * )  *  3 1  4 1 6  9  3  F  »0  .5oOU,"Jb  00 

AMD 

C. 2559665b  0i-U.B785684b-02-').26*4627C  02 
RH 


o.c 

6.5000 

1.0 JOO 

l.50o  » 

2 . SoOO 

3 . 01  V 

i.'A  ^ 

4.  5 ''Ob 

5.0000 

lO.OCuO 

5.5000 

6 • 009 u 

6. SO  )  ’ 

7.  <■‘00. 

7.  S 

p.oor: 

8  .SCO  ) 

9.  DO'lO 

9.S.  0'' 

ZH 

O.C 

0.  »J 

c.o 

0..-' 

O.C 

0.0 

c.o 

0.0 

0.0 

c.{ 

o.c 

O.C 

o.u 

(j  .  u 

0.0 

0.1 

9*  J 

C.u 

0 . 9 

C.o 

o.f 

ELECTRIC  HELD  ANO  i,AIN  OF  U  Ol'it  Jk«*.  \T  l\U  LAM30A  =  <..26^  *1 


t# 

REAL ( EG ) 

IMAM  E9) 

GA  [  Nb> 

RbAL (  L/ ) 

I MAb( htf ) 

GAIN/ 

c.o 

-C . U04b 

0  1 

0.  ) 

0.L699C  01 

0. 1 304F 

91 

O.C 

9. 1499b 

r  l 

10. 0 

-0. L243t 

01 

9.3 

0.  15461-  01 

0.  1134b 

ul 

3.0 

J. 1287b 

;  i 

2C.U 

—  0 • 1 08  uE 

)l 

0.0 

<'.11'-51  01 

3.68C3E 

91 

3.0 

" .4628b 

37.0 

-0. 8*3 5 8t 

ro 

0.  ) 

1, .  7  32  4F  J  : 

0.  74  77b- 

31 

3.9 

0.6501E- 

.  2 

41i*  0 

-C.6218E 

cu 

0.0 

<-.J8(>7t  1C 

-0.5131b 

>  L 

3.  ' 

'./842b 

!  * 

50.0 

-0. 4158t 

00 

0.0 

C.1729F  00 

-#:#lr36C 

0) 

3.  9 

j.  U  73 E 

C  1 

6r  •  0 

-0.2558E 

00 

0.0 

C.6S42F-91 

-0.  1391b 

01 

3  •  0 

9 . 1 9  3  5  F 

:  i 

7  J  •  0 

-0.1416b 

00 

0.0 

f.2005fe— J1 

-C.  1607  b 

Cl 

3.0 

0.258/b 

1 1 

8:  .  0 

-0 .61 32b- 

91 

0.0 

'■>.  1822E-02 

-I/.1716F 

Cl 

0.0 

C./944E 

C  l 

90 . 0 

-0. 1056b- 

•06 

>. 9 

U. tll«c-l3 

-o. 1748C 

91 

0.3 

~.3"r>7b 

r  l 

LOw.O 

C  .6182E- 

■01 

3.  9 

0. 3«22b-02 

-0. 1716F 

91 

0.0 

C.2944F 

;  l 

113.0 

0.1416b 

e  j 

0.3 

0.2C0Sb-'.U 

-0. 1607b 

Dl 

0.0 

0  •  2  5  8  ?  F 

C  1 

1.20.0 

0.2558b 

00 

0.3 

0«  6542F-01 

-0.1391b 

01 

0.0 

9 . 1 934  E 

01 

130.0 

0.4158E 

03 

3.0 

0.1729F  jC 

-0. 1036b 

01 

0.^ 

0.1  . 73^ 

01 

140.0 

0.621 BE 

oO 

0.0 

0 . 386 7t  00 

-C.5331C 

oc 

*9. 0 

v .2842b 

i  ' 

150.0 

0.8658b 

CO 

0.3 

c. 7324b  00 

0. 7477b- 

■01 

3.0 

j.54^ir- 

160.0 

0.  ioBOE 

01 

3.0 

*.LlG5b  31 

*j  %  68  03  f 

ct 

0.0 

C*  .4  u?  8E 

170.0 

C.1243E 

01 

0.  9 

0.1546F  01 

0.1134t 

ul 

3  •  u 

r  .1287b 

r  i 

18C.0 

0.1304E 

01 

0.0 

0.1699?  01 

0. 1304 F 

01 

3.0 

M  1699E 

c  1 

ELECTRIC  FIELD  AND 

GAIN  HF  L  IGfJCdKBUiT  F'JS  LAMrtDA  = 

?—•  • 

8736T-  2 

S 

P  CAL l fc£ ) 

IMAM  19) 

(i  A  I  Nc 

RbAL ( E/I 

I  MAM  b'/  ) 

n  ain/ 

0.0 

-0.2266E 

J1 

u*  ) 

C ♦ 6 1 36b  j 1 

C.  2266b 

91 

0  .  u 

* ' .  5  1 3  6  E 

r  l 

10.0 

-C.2L76E 

01 

0.3 

C.4/36F  01 

0  .  i  1 8C  t 

01 

3.0 

0 .475  2 1 

i:  i 

2  J  .  0 

-0.1929b 

01 

0,0 

1)  ♦  3  7  2  2  b  91 

0.  1943F 

Ul 

3.0 

0. 3776b 

r  1 

30.0 

-0. 1534E 

01 

0.3 

0.2510b  'J  l 

J. 1614F 

t  1 

3.0 

0.2637E 

:  i 

4y  .  0 

-0.1210b 

)l 

9.0 

0.1465b  01 

0.1263b 

01 

3.0 

) .  1695F 

*  1 

6).0 

-C. 8632b 

00 

0.  ) 

0  •  7452t  00 

0.9445b 

or 

0.0 

0.8O21F 

*1 

60.0 

-0  •  6  722b 

00 

0.  9 

0.3275E  00 

0.6923b 

Ol 

0.0 

0 .47930 

V 

7  3*.  0 

-0.3416b 

00 

0.) 

C.1167E  00 

0.5167F 

Ou 

O.r 

0  •  2  6  7  r  fc 

c 

80.  C 

-0 . 15S2b 

0) 

0.3 

0. 250  lb-0  1 

0.4154F 

00 

O.C 

0. 1726b 

'  r 

V 

9C.0 

-G.2768b^06 

0.0 

0. 76611-13 

•i.3325b 

oc 

0.0 

M  14630 

* 

100.0 

0.1582b 

CO 

0.0 

0.25MF-G1 

0.4154c 

or 

0  •  G 

'.1/250 

110.0 

0. 3416E 

30 

0.9 

C. 1167b  00 

0.5167F 

oc 

0.0 

0.2O7CF 

V 

123.0 

0.5722E 

00 

C.O 

0.3275F  00 

0.6923b 

oc 

O.U 

0.4793b 

rt 

130.0 

0.8632* 

00 

o.y 

C.7462F  00 

0.9445E 

oc 

0.0 

0.8921b 

>  ** 

140.0 

C.121CF 

01 

0.3 

0.14o5r  01 

0. 1263b 

01 

0.0 

0.  1695F 

0  l 

150.0 

0.  1584C 

01 

0.0 

C  »  25 1  OF  01 

0. 1614F 

01 

0.0 

L.?607f 

r  l 

16C  .  C 

0.1929E 

91 

0.9 

0.3722b  01 

0. 1943 F 

01 

0.0 

0.3776E 

r  i 

170.0 

0.2176E 

31 

0.3 

0.4?  36b  01 

0.21 8ftfc 

01 

0.0 

U.4752F 

f  l 

180.0 

0. 2266b 

01 

0.0 

G.5136E  31 

0.2266F 

01 

0.0 

0.5136F 

Cl 
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fcltCTKlC  HtlU  AND  t.AIM  (U  1  1GFNC JRkFnT  FU«  LAMBOA  -C-.?6'obf  P2 


p 

RtAL(fc6) 

I  MAG ( L  0 ) 

6AING 

KCALIh'/) 

l a Ab( b  ) 

•  J 

C. 1836b 

31 

0.0 

0.3371b 

01 

-0.1836b 

Cl 

0.0 

0.337 IF  •.  1 

i  o  .  v 

C.  150 Jt 

01 

0.0 

0. 2260b 

Ul 

-0. I 7  2  3  r 

1*1 

0.0 

o.poftop  :i 

2  '  >  •  V 

0.6433b 

GO 

3.0 

v,.4l  39b 

JO 

-0. 14196 

01 

0.0 

0.2-1  At  .  i 

3J.G 

-C.4076E 

00 

0.0 

0 . 1 66  i  r 

30 

-0.L011C 

01 

0.0 

0.1  >?16  '  1 

4C  . 0 

-0.1289b 

01 

0.0 

0. 160  IE 

01 

-0.5953b 

oc 

Oc  0 

0  •  3*>44F  r- 

5..1 

-0. 1  7666 

Jl 

0.  J 

0. 3120L 

01 

-0.2446b 

00 

D .  9 

r.$9AAE-<  1 

6^.0 

-0.1789b 

01 

0.0 

0.3169b 

Jl 

0.1034b- 

■01 

3.0 

•  KhAt  i 

7).l 

-C. 1403k 

31 

3.0 

19081 

Ul 

C.  17  16b 

Ot 

0.0 

5.?)4&f-:  < 

H".P 

-C. 76426 

VO 

J.) 

C  ♦  584 'T 

JO 

0  •  2  8  6  8  b 

oc 

3.0 

'  .6S0  31-C  1 

90,0 

-i.  .1412E- 

05 

0.0 

0.1992b- 

-ll 

0.28296 

•A 

0.0 

G.S  'Obt-C  1 

IOC .  r 

0 .7642F 

JO 

0.3 

0.5840b 

00 

0.2568T 

00 

D.C* 

0 .6r>93t-r  1 

1LC.C 

0. 1403F 

01 

0.0 

0. 1963F 

01 

0.1716b 

91 

0.0 

0.?94bF-; l 

12.  .  > 

C  •  l  78  )fc 

01 

J  .  J 

0.3101 

n 

0. 1  u  3  4  f  - 

■01 

O.u 

G.  K.VIF-'*  ' 

13^  .  „ 

C • 1 76ofc 

)1 

0.  ) 

0.  31  2  )b 

oi 

-0.2446b 

no 

0.0 

0. 303AC-: L 

*  u 

C. I20)t 

ul 

0  .  ) 

u  •  L  6o  1 F 

01 

-0 • 5953b 

0»‘ 

3.0 

J.33AAF  •* 

15). 0 

G. 40766 

00 

0.0 

C. 1661b 

ce 

-0. 10 1 lb 

»:i 

0.0 

0. 102  IE  '1 

16  ) ,  0 

-0.64336 

0  3 

0.0 

C. 4139b 

u( 

-0.1419b 

91 

0.0 

0.2O14F  ;  l 

17*: 

-0. 1503b 

Jl 

0.0 

0  o  226GF 

01 

-0.1723F 

01 

0.0 

0.29(>9t  1 

18". 0 

-C.  1536b 

Jl 

0.  ) 

0.33716 

)i 

-0.1836b 

Al 

3. 
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VI.  SCATTERING  CROSS  SECTIONS 

Program  #5  which  calculates  and  plots  o/A2  for  an  axially  incident 
plane  wave  accepts  punched  card  data  according  to 
READ  (1,10)  NP,  NT,  NS,  JM,  BK 

10  FORMAT  (413,  E14.7) 

READ  (1,11) (AMD (I),  1=1,  JM) 

11  FORMAT  (5E14.7) 

READ  (1,15)  (RH(I) ,  1=1,  NP) 

READ  (1,15) (ZH(I),  1=1,  NP) 

15  FORMAT  (10F8.4) 

READ  (1,50)  (L(T.) ,  1  =  1,  JM) 

50  FORMAT  (2013) 

The  variables  NP,  NT,  NS,  JM,  BK,  RH  and  ZH  are  the  same  as  those  in  program 
it 4.  The  variable  AMD  appears  under  the  heading  "eigenvalues  of  the  matrix  B" 
in  the  output  of  program  it 2.  The  L(I)  eigencurrent  is  the  I  eigencurrent 

to  be  considered  in  the  modal  expansion  of  the  scattered  field.  The  variable 
L(I)  is  necessary  because  it  is  desirable  to  perform  the  modal  expansion  by 
adding  eigencurrents  in  order  of  increasing  |  A  | .  (Program  // 2  has  ordered 
the  eigencurrents  in  order  of  increasing  A.)  For  instance,  if  A^  corresponds 
to  the  smallest  |a|  and  X<-  corresponds  to  the  next  smallest  |a|,  then  L(l)  =  6 
and  L(2)  =  5.  The  impedance  matrix  and  eigencurrents  are  read  from  records  1 
and  2  of  direct  access  data  set  6  according  to 

REWIND  6 

READ (6) (Z (I),  I  o  1,  NZ) 

READ (6) (FI (I) ,  1=1,  NZ1) 

where 

NZ  =  NM2*NM2 
NZ1  =  HM2*JM 

and  NM2  is  either  NP-i  or  NP-3  depending  on  whether  or  not  the  generating 
curve  C  closes  upon  itself. 


32 


Che  subroutine  PLANE  is  the  same  as  in  program  // 4.  Much  of  the  logic 
before  statement  7  in  the  main  program  is  devoted  to  preparing  the  input 
data  for  plane  and  is  thus  the  same  as  in  program  // 4* 

DO  loop  30  modifies  the  impedance  matrix  according  to  (2-47)  of  [1]. 
DO  loop  19  obtains  the  matrix  X  appearing  in  (2-15)  of  [1]*  Statement  6 
inverts  the  impedance  matrix  to  obtain  the  admittance  matrix  Y. 

DO  loop  85  calculates  the  scattered  field  by  inserting 


rt-i 


LJ*  J 


=  [I]  =  [Y]  [rJ 


(15) 


into  (5)  for  n=l.  Equation  (15)  is  possible  because  plane  wave  excitation 
and  measurement  coefficients  have  the  same  form*  Notice  that  (1-48)  of 
[1]  is  written  using  the  eigencurrents  as  a  basis  while  (15)  is  written  with 
the  expansion  functions  (2-42)  in  mind.  The  excitation  of  a  u^  polarized 
plane  wave  incident  in  the  plane  $  =  0  gives  rise  to  expansion  functions 
(2-42)  of  [1]  while  the  excitation  of  a  polarized  plane  wave  in  the 
plane  $  =*  0  gives  rise  to  expansion  functions  (2-43)  of  [1]*  In  DO  loop  85, 

K5=l  obtains  plane  wave  excitation  from  the  direction  (0=tt,(J)=O)  and  K5=2  that 
from  the  direction  (8=0,4>=0).  DO  loop  82  stores  the  column  matrix  [I]  of 
(15)  in  E3.  DO  loop  103  stores  E0,  E^,  (a/A2)^,  and  (c/A2)^  in  E(J),  E(J4NT), 
SIG(J)>  and  SIG(J+NT)  where  J  indicates  the  J^1  measurement  polar  angle  0.  Here, 
is  the  0  component  of  the  far  zone  scattered  field  in  the  plane  $  =  0  and 

IT 

Ea  is  the  <p  component  of  the  scattered  field  in  the  plane  <}>  =  -r.  The  phase  of 

t  —  i  ler  ^ 

the  scattered  field  is  normalized  to  -je  J  *  fhe  magnitudes  |Eg|  and  |E^|  are 


such  that 


(oA2)e=|Ee|2 
<«A2)d-|E  J2 


(16) 


The  constant  Cl  used  to  normalize  the  scattered  field  is  given  by 
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The  factor  w2u2/ (4xX2)  appears  in  (1-51)  of  [1],  The  logic  between  state¬ 
ments  119  and  109  finds  SCL(K5)  so  that  the  maximum  value  of  SIG(J)*SCL(K5) 
for  J  =  1,2,...2*NT  is  always  between  3.2  and  3.0.  DO  loop  106  prepares  hori¬ 
zontal  and  vertical  coordinates  E5  and  E6  suitable  for  plotting  o/ A2. 

DO  loop  48  multiplies  the  eigencurrents  by  p  to  retrieve  the  coefficients 
[I]  of  the  expansion  functions  (2-42)  of  [1].  The  matrix  [I]  is  stored  in  FI 
with  subscripts  (J-1)*NM2+1  to  J*NM2  denoting  the  eigencurrent .  The  admit¬ 

tance  matrix  [Y]  appearing  in  (15)  will  be  expressed  in  terms  of  [I]  whose 
columns  are  the  expansion  coefficients  of  the  eigencurrents.  For  the  moment 3 
assume  that  there  are  NM2  eigencurrents  so  that  [I]  is  a  square  matrix.  Equa¬ 
tion  (2-18)  of  [1]  can  be  multiplied  by  [I]  to  yield 

[I]fZ][I]  *  [I] [X] [I]  (j  +  l/A)  (18) 

Because  of  the  orthogonality  relationships  (C-24) ,  of  [1],  the  right  side 
of  (18)  is  a  diagonal  matrix.  Equation  (18)  can  be  inverted  to  obtain 

[X]_1[Y] [I]”1  - - - -  (19) 

tiKxHi]  (j  +  1A) 

which  leads  to 

m  =  [i] - - -  [i]  (20) 

(inxni]  (j  +  1A) 

Equation  (15)  is  introduced  into  (5)  with  the  result 
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where  [Y ]  is  given  by  (20).  Equation  (21)  is  a  specialization  of  (1-37) 
of  [ 1 3 *  Alternatively,  (21)  could  have  been  obtained  by  finding  the  matrix 


[IC] 

il*') 


associated  with  J  of  (1-30)  of  [1]  and  inserting  this  matrix  into  (5). 


Evidently,  the  use  of  JM<NM2  terms  of  the  sum  (1-37)  of  [1]  is  equivalent  to 
using  only  JM  columns  of  [I]  in  (20).  DO  loop  21  stores  the  diagonal  matrix 

[1  ]  ['X ]  [ I  ]  ( j  +  l/\j  aPPearin8  in  (20)  in  E3.  The  inner  DO  loop  27  stores 
1 


[I] 


[I][X][I]  (j  +  1/A) 


in  T3. 


The  index  K  of  DO  loop  29  indicates  that  K  columns  of  [I]  are  being 
used.  DO  loop  30  adds  the  contribution  of  the  L(K)C^  column  of  [I]  to  the 
admittance  matrix  [Y],  In  DO  loop  65,  K4  =  1  obtains  the  plana  wave  incident 
from  (0=n,  4>=0)  and  K4=2  that  from  (0=0,  <J>— 0) .  Inner  DO  loop  45  stores 

1^1  [R^  -j ]  in  E3.  In  DO  loop  66,  K5=l  obtains  E^  and  K5=2  obtains  E^. 

DO  loop  44  calculates  E„,  E.,  (o/A2)rt,  and  (a/ A2),  and  stores  them  in  E  and 

b  9  fc;  9 

SIG.  The  quantity  o/A2computed  using  the  modal  approximation  to  [Y]  is 
plotted  as  symbols  (X  for  0  polarization  and  0  for  <j>  polarization)  and  the 
quantity  of  A2  computed  using  [Y]  obtained  by  inverting  the  impedance  matrix 
is  plotted  as  a  solid  curve  in  the  latter  part  of  DO  loop  66.  The  logic 
following  statement  66  is  devoted  to  drawing  the  axes  for  the  plot  and  print- 

The  scattered  field  components 

o  V  o  V 

E, 


ing  the  quantities  E^,  E^,  (a/A2)^  and  (o/A2)^. 

E,  that  are  printed  lack  the  phase  factor  -je  and  are  normalized  accord- 


'6  ’  "6 
ing  to  (16) . 


Minimum  allocations  are  given  by 

COMPLEX  Z(NM2*NM2) ,  Y(NM2*NM2) 3  VR(2*NT*NM2) ,  E3(NM2),  E(NT*2),  T3(NM2*JM) 
DIMENSION  AMD(JM) ,  RH(NP) ,  ZH(NP) ,  FI (NM2*JM) , 

DH(MP-l) ,  TH(NT) ,  X(NM2*NM2) ,  SIG(NT*2) , 

L (JM) ,  SN(NT) ,  CS(NT) ,  E9 (NT) ,  E10(NT), 

E5(NT*8) ,  E6(NT*8) 


COMMON  RS(NP-l),  ZS.(NP-l)  ,  SV(NP-l),  T(NM2*2),  TR(NM2*2) 

DIMENSION  B J ( 3*  (NP-1)) 

DIMENSION  LR(NM2) 

Here,  NP  is  its  value  after  execution  of  statement  90  in  the  main  program. 
Note  that  BJ  appears  in  PLANE  and  LR  in  LINEQ. 
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Listing  of  Program  it 5 

//  (0034, EE, 4, 2, ,7) , 'MAUTZ , JOE • ,MSGLEVEL=1 

II  EXEC  FORTGClGtPARM.f-riRT^MAP' 

//H1RT.SYSIN  DD  * 

SUBROUTINE  PLANE ( VVR , THR , NT ) 

COMPLEX  VVR ( 1 ) , A5  *  A6 » U 

COMMON  U,RS(40) , Z  S (40 ) » SV ( 40) ,CV (40) , 8K»NP,NN,T (80) , TR ( 80 ) 
0 1  MENS  I ON  R.M  126) »THR ( 1 ) ,FK( 20) 

KG=NP-1 

NM=KG/2-l 

M2=NN+2 

A5  =  2.*3.141593*IJ**(NN+1 ) 

NV=NM*4 
FK ( 1 ) *1  • 

DO  153  J=1 *M2 
J1=J+1 

FK( J1)=FK( J)*J 
153  CONTINUE 

DO  156  L=1,NT 
Ll= ( L-l ) *NV 
CS=COS ( THR ( L ) ) 

SN=S I M ( THR ( L ) ) 

BCS=BK*CS 
DO  302  J= 1 , KG 
X=RS 1 J ) ^BK^SN 
J1  =  J 
I  1=NN 

IF(I1)  303,304,303 

304  11=11+1 
J 1 = J 1+KG 

303  DO  305  JJ=I  1,M2 

I F ( X-l . E-5 )  1,1,2 

1  IF(JJ-l)  3,3,4 

3  BJ ( J1 )  =  1  • 

GO  TO  306 

4  BJ ( J1 ) =0. 

GO  TO  306 

2  RH=X/2. 

RH2=RH*RH 
RH3=RH**t  (  J  J— 1 ) 

BJIJ1)=RH3/FK(JJ) 

SS=BJ ( J 1 ) 

ft  SST=SS*1.E~7 
DO  155  K= 1 , 20 
SS=-SS*RH2/K/ (  K+  J  J-l  > 

BJ)  J11=BJ(  JD+SS 
I F ( ABS ( SS ) -SST )  306,306,155 
156  CONTINUE 
STOP  155 

306  J1=J1+KG 

305  CONTINUE 
302  CONTINUE 

IF(NN)  307.308,307 

308  DO  309  J= 1 , KG 
J]=J+2*KG 

BJ ( J ) =-B J ( J 1 ) 

309  CONTINUE 

307  DO  300  J=1 ,NM 
J 1=J+L 1 

J2= Jl+NM 
J3= J2+NM 
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J4= J3+NM 
VVR ( J 1 ) =0 . 

VVRl J2)=0. 

VVR ( J3 ) =0. 

VVR ( J4 ) =0. 
on  301  1-1,4 
I 1=2*( J-l J+I 
1 4=4;*  (  J-l  )  +  l 
I2»U+KG 
1 3= 1 2+KG 

A6  = (COS ( 7  S  ( I  1 )*RCS)+U*SIN(ZS(I1 )  *BCS ) )*A5 
R.)1=(BJ(  I  3 )  +BJ  ( I  1 )  )  #  .  5 
BJ2=(B.U  1 3 ) -P J ( II )  )#.5 

VVR( J1)=VVR( J1)+AA*(CS*SV( I 1 )  *B  J2+SN*CV  (  II  )*BJ(  I2)*U)*T( In) 
VVR  (  J2 )  =VVR  (  J2 )  +A6*CS*BJ1*TR  (14) 

VVR(  J3)=VVP.(  J3)  +  A6*SV(  I1)*BJ1*T{  14) 

VVR ( J4 ) =VVR ( J4 ) +A6*BJ2*TR (14) 

301  CONTINUE 
300  CONTINUE 
136  CONTINUE 
RETURN 
END 

SUBROUTINE  L I NEO ( LL ,C ) 

COMPLEX  Cli) ,STOR,STO,ST,S 
DIMENSION  LR(5B) 

DO  20  1  =  1, LL 
LRU  )  =  I 
20  CONTINUE 
M1=0 

DO  18  M= 1 , LL 
K=M 

DO  2  I =M ,  L  L 
K1=M1+ l 
K?=M1+K 

I F ( CARS ( C ( K 1 )) -CABS  I C ( K2  )  )  )  2,2,6 

6  K=  I 

2  CONTINUE 
LS=LR(M) 

LR(M)=LR(K) 

LR ( K ) =LS 

K2=M1+K 

ST0R=C(K2) 

J1  =  0 

DO  7  J  =  1 , LL 
K1=J1+K 
K2=J 1+M 
ST(1=C(K1  > 

C(K1)=C(K?) 

C(K2)~ST0/ST0R 

J1=J1+LL 

7  CONTINUE 

K*  ,IM  ■  It 

l  — n  i  Ti*! 

C(K1)=1./ST0R 
DO  11  1  =  1, LL 
IF(I-M)  12,11,12 
12  K1=M1+ I 
ST=C(K1) 

C ( K1 ) =0. 

J1  =  0 

DO  10  J=1 ,LL 
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k  ]  =  Jl  +  I 

K  ?  =■  J  1  +M 

r  (fvl)=C(Kl)-f  ( K  2 ) '  ST 
J1=J1+LL 
in  COMT  I NOF 
ll  CONTINUf 

mi=mi+li 

18  COM  I  NOE 
J1  =  0 

no  9  .!=i, ll 
mj-LRun  14,8,14 

14  L  8  J  =  L  K  (  J  ) 

J2=(LRJ-l)*Ll 
21  OH  13  1=1  *  LL 
K2=J2+I 
K  1  =  J  1  +  1 
S=C(K?) 

C ( K2 ) =C 1 K  1 ) 

C  (  K  1 )  =  S 
13  CONTINUE 

LR« J)=LR(LRJ) 

LR!LRJ)=LRJ 
I F ( J-LR ( J ) )  14,8,14 

ft  J 1= J 1+LL 
9  CONTINUE 
RETURN 
END 

COM  PI  EX  U,01,Z  ( 1444)  ,Y< 1444)  ,VR(  554ft)  ,E3(38)  ,  6  (  146  )  ,  T3  ( 1444 ) 

0 1  MENS  1 ON  AMD  1 3ft 1 , RH ( 4 1 ) , 2 H ( 4 1 ) , F I ( 1444 ) , DH ( 40 ) ♦ TH ( 73 ) , X ( 1444 ) 
DIMENSION  SIGl 14  6)  ,L(3ft) , AREA (400) ,SN( 73) ,CS ( 73 ) ,E9 ( 73 ) ,E 10 ( 73 1 
DIMENSION  XP(2) *YP(2) »E5(584),E6(584) ,SCL(2) 

EQUIVALENCE  ( Z (  1  )  , Y ( 1 ) ) 

COMMON  U»RS ( 40 ) , ZS ( 40 ) » SV(40) »CV(40) » BK»NP»NN,T ( 80) ,TR ( 80 ) 

CALL  PLOTS! AREA, 400) 

READ! 1 , 10)  NP,NT.NS, JM,BK 

10  FORMAT (4 13, El 4. 7) 

READ! 1, 11) (AMD! I ) , 1=1, JM) 

11  FORMAT! 5E 14. 7) 

READ! 1, 15 i !RH( I ) , 1=1 ,NP) 

READ! 1,15) ( Z  H ( I ) , I  =  1 »NP ) 

15  FORMAT! 10F8. 4) 

READ! 1,50) (L( I) , 1=1, JM) 

50  FORMAT (2013) 

WRITE (3,9) 

R  FORMAT! '1  NP  NT  NS  JM  BK*  ) 

WR I T E ( 3 , 1 2 )  NP ,NT ,NS * JM,BK 

12  FORMAT ( IX, 413, E 14. 7) 

WR I T  F ( 3 , 1 3 ) (AMD! I ) ,1=1 ,JM) 

13  FORMAT! 1 OAMD ' / ( 1X,5E14.75 ) 

WR (TE (3,16) (RH! I ) , 1=1 ,NP) 

16  FORMAT! 1 ORH ' / ( lX,10Ffi.4) ) 

WRITE! 3, 127? ! ZH ( I )  ,1=1  ,NP) 

127  FORMAT! 'OZH'/( IX, lOFft.4) ) 

WR I T  E ( 3 , 5 1 ) ( L ( I ) ,1=1 , JM) 

81  FORMAT! '0L'/( 1X.20I3) ) 

P]=3, 141593 
E  T  A  =  376. 730 
O=(0.,l.) 

Cl  =  RK’-"ftK#E1  A/(4.*S0RT(Pl**3)  ) 

KL  =  1 
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I F ( (RH( l)-RH(NP) ) .NE.O. .OR. (ZHi 1 )-ZH{NP) ) .NE.O. )  GO  TO  90 
KL  =  0 

RH ( NP  + 1 ) =RH ( 2 ) 

ZHtNP+l )=7H(?) 

RH ( NP+? ) =RH ( 3 ) 

ZH(NP+2)=ZH(3) 

NP=NP+2 

90  NM2=NP-3 
NZ=NM2#NM? 

NZ 1=NM2*JM 
RFWINO  6 

READ ( 6 ) (Z (!)  ,  1=1  ,NZ  ) 

READ (6) (FI ( I ) ,1=1, NZ1 ) 

NM=NM2/2 
NM4=MM*4 
NT2=NT*2 
NT3=NT-NS 
NT4=NS+1 
DO  40  I=?,NP 
12=1-1 

RR 1=RH ( I ) -RH ( 12) 

RR2=ZH( I )-ZH( I?) 

DH (  I2)=S0RT(RR1*P.R1+RR2*RR2) 

RSt I2)=.5*(RH( I )+RH(I2)  ) 

ZS(  1 2 )  = . 5* ( ZH (  I  )  +7.H  (  12)  ) 

SVU2)=RR1/DHU2) 

CV(  I2)=RP.2/DH(  1 2 ) 

40  CONTINUE 

DO  41  J=1,NM 
J2=2* { J-l ) + 1 
J3= J2+1 
J4= J3+1 
J5=J4+1 
J6=4* ( J-l ) + 1 
J7= J6+1 
JB  =  J74-1 
J9= J8+1 

DEL1=DH ( J2 ) +DH ( J3 ) 

DEL2=DH ( J4 ) +DH ( J5 ) 

T  (J6)=DH(J2)  ^DH  (»)2)/?,/DELl 
T(  J7)=DH(  J3)*(DHU?)+DH(  J3)/2.  )/DELl 
T ( J8 ) =DH ( J4 ) v ( DH ( J5 ) +DH < J4 ) /2  .  ) /DEL2 
T( J9)=DH( J5)*DH( J5) /2./DEL2 

41  CONTINUE 

DO  91  J= 1 i NM4 
TR( J)=T( J) 

91  CONTINUE 

I F I KL  » ED, 0 )  GO  TO  95 
I F  <  RH ( 1 ) ) 93  *  94 , 93 

93  DEL1=DH 1 1 ) +DH ( 2 ) 

TR( 1)=DH( 1)*( l.+(DH(2)+DH( l)/2. )/DELl) 

TR(2)=DH(2)*t l.+DH(2)/2,/DELl ) 

94  IFIRH(MP) )96,95,96 
96  Jl= ( NM-1 ) *4+3 

J2  =  .ll  +  1 

DEL2=DH(NP-2)  +DH (NP-1 ) 

TR { J 1 ) =DH( NP~2 ) * ( 1 ,+DH (NP-2)/2./DEL2) 

TR ( J2 ) =DH ( NP-1 )* { 1 .  +  ( DH { NP-2 ) +DH (NP~l)/2.)/DEL2) 

95  DEL=PI / ( NT-1 ) 

DO  43  J=1,NT 
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1H|  j)  =  (  j-1  UDEL 
SN( J)=SIN(TM(J) ) 

CS(  J)=C0M  IHU)  I 
CONTI  NIN¬ 
ON  =1 

7  CALL  PLANE(VR,TH,NT) 
0EL1=1«0./Pl 
00  12P  J=l,NT 
TH(  J)=TH(  J)*[)6L1 
128  CONTINUE 
XP! 1 )~2r 
XP ( ? ) -8  . 

Y  P ( l ) =5 . 

YP ( 2 ) =5 . 

01  =  .5<‘U 
00  80  4=1, NM 
Jl=( J-l)*NM?+NM 
J4=( J*MM-1 )*NM2 
00  81  I  =  1 , NM 
J2= J 1 + I 
J8=J^+I 

Z IJ2)=U1*Z( J2) 

7 ( J3I=-U1*Z ( J3) 

J5=J?-NM 
J6= J3+NM 
Z(J5)=.5*Zi J5) 

Z( J6)=.5*ZI J6) 

PI  CONTINUE 
80  CONTINUE 

00  19  J~1  ,NM2 
J2=U-1)*NM2 
DO  20  1=1, J 
J3= J2+ I 

J4=(  J-1  KMM2+J 
Ul=.b*(Z( J3I+ZI J4)l 
X  (  J  3  )  =  A  I M  AG  (  U 1  I 
X( J4)=X( J3) 

20  CONTINUE 
19  CONTINUE 
ft  CAIL  L  I NEO  ( NM2  ,Z  ) 

DO  85  K5= 1  ,  2 

J5=  (  2-K5  )  "7  ( NT-1 )  *NM4 

00  R 2  1  =  1, NM2 

E3II)=0. 

DO  83  KK«1»NM2 
J3=I+(KK-1)*NM2 

J4= J5+KK 

E  T { I ) =E3 ( I )+Z (43) *VR ( J4 ) 
83  CONTINUE 
R?  CONTIMUF 
S2  =  0. 

00  103  Kft=l,2 
J5=( Kfi-1 5*NM2 
J8=(Kft-l  I-XNT 
00  84  J= l .NT 
0  7  =  J  +  J  8 
111  =  0. 

Jl=( J-I )#NM4fJ5 
DO  88  1=1, MM2 
J2=J1+! 
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U1=U1*VR(  J?)*F.3(  I  ) 

88  CONTINUE 
F  (  J7  )=C  1*1)1 
Sl=CAB$fE( J7  )  ) 

S I G { J7)=S1*S1 

lr  (  S I G {  J7)  .GT.S?)  S2-S1GU7) 

84  CONTINUE 
103  CONTINUE 

WR  I  T  h  {  3 , 5  ) 

5  FORMAT (  1 OSCATT  F  &  En  FIELD  AND  SCATTERING  CROSS  SFCT ION/WAVEL ENGT H  S 
10UARFD 1 ) 

Sl=( 2-K5)*180. 

WRITE(3,124)  SI 

124  FORMAT ( 1  BY  MATRIX  INVERSION,  INCIDENCE  FROM  0=',F4.0,'t  0=0') 

2  WRITE (3, 118) 

11R  FORMAT (  '  +  *  1 36 Y , 1  - 1  ,7X,  */») 

WRITE (3, 39) 

39  FORMAT  (  1 0  0  R  t  A  L  ( b  0 )  IMG(EO)  SO/ILAM)**?  RFAL(FO) 

llMAG(EO)  SO/ ( LAM)##? 1 ) 

WRIT E(3,117) 

117  FORMAT (  *  +  -* « 1 IX,  ■-  1  ,  11X, - « , 1 8X , • / • , 1 1 X ,  *  /  /') 

DO  37  J=ltNT,NS 
J 1 = J+NT 

WRITE  (3, 38)  TH(J)  ,E(  J)  ,SIG(J)  tE(  Jl)  tSIGUl  ) 

38  FORMAT ( IX, F6. 1 , 6E 12*4 ) 

37  CONTINUE 
119  J 1= 10+AL0G10 ( S2 ) 

S3=.l**( Jl-10) 

$4=S2*S3 

I  F  C  S4-1.5)  no, no, 111 

110  $Cl(K5)*2.*S3 
GO  TO  112 

111  I F { S4-3  .  )  113,113,114 

113  SCL  \  K5  )  -  f,3 
GO  TO  11? 

114  I F ( $4-6 ♦ )  115,116,116 

115  SCL ( K5 ) - • 5*S3 
GO  TO  112 

116  SCL(K5)»«2*S#* 

112  S5= 1 ♦ /SCL ( K 5 ) 

WRITE (3, 109)  S5 

109  FORMAT ( 1 OOME  INCH  CORRESPONDS  TO  S I GMA/ ( L AMROA ) £#2= 1 ,E11*4) 

DO  106  K6= 1 , 2 
J 1=  (  K6-1 )  *’  NT 
J3=NT2*< ( K5- l ) *2+K6— 1 ) 

DO  107  J= 1 ,NT 
J  7  =  J 1  +  J 
J  8  =  J  +  J  3 
J9= J3+NT2-J+ l 
S2=S1G( J7)*5CL(K5) 

S1=SU( J)*S2 
E5U8)=5.  +  S1 
E5( J9)=D.-S1 
E6(  J8) -5*+CS( J)#S2 
E6U9)=E6(  JR) 

107  CONTINUE 
106  CONTINUE 

85  CONTINUE 
DU  48 

J1=C J-1)*NM2 
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>IM  ]8  1  =  1,  NM 

.M-Jl  +  I 

.1  *=•  J.>  +  M%1 
J4~?‘  1  +  1 

E  I  (  22 1  =l  I  ( ,1 7 1  .-fin  ( .14  ) 

H  C  J31  =H  I  I3)«PH(.|4) 

18  CONTINUE 
4«  cpntinui- 
nu  21  J= 1 , JW 
Jl=» J-l >  -Nf->? 

Sl=n. 

IM  22  1=1 ,NM2 

S2=0. 

J4=  t  I -1 )vNH2 
lid  23  K=1,NM? 

J3=J1+K 
J2= J4+K 

52  =  S2  +  X  ( .12  )  ;:F  I  (,I3! 

23  CUNT  I  Ml  IF 

.12  =  J 1  +  I 

Si=Sl+S2“F I ( J2) 

2?  CONTINU¬ 
ES!  J)  =  1./S1/(II+1./AHD(J)  ) 
DP  27  1=1 ,NM2 
J2= J1 + 1 

T  3  (  J  2  )  =  r  I  ( .12  )  >;63  ( -1 ) 

27  CONTINUE 
21  CONTINUE 

•IP  32  J=1,N? 

Y  ( J )  sfl, 

32  CONTINUE 

1)0  29  K  =  l,.IM 
52  J3= ( L ( K )  —  1 )*NM2 
OP  30  J  =  1  , MM  2 

.ii= (  j-h*mm2 
J5= J3+J 
UO  31  l  =  l,.l 
J2= J 1+ I 
J4= J3+ 1 

Y(  J2)=Y(.I2)-fFI  ( J5)*T3(.j4) 
26=  ( 1-1 )#MH2+J 

Y  ( ,16  )  =Y  ( .12  ) 

31  CONTINUF 

30  CONTINUE 
no  65  K4=l ,2 

J8  = ( 2-K4 ) *  1  NT -1 1 =NW4 
HP  45  1=1,NM2 
E3(!)=0. 

Ilf)  46  KK  =  l,NM2 
J3=I  +IKK-1  )>!NM2 
J4  =  J.3+KK 

F3 1  I )  =E3 ! I  5  tY I J3) YVR I J4 ) 
46  CONTINUE 
45  CONTINUE 

00  66  K5=l .2 

K3=  I  {  K 4  —  1  )  ><2+K5- 1  1  >'NT2 

,17=  ( K  5”  1  )*NT 

K2=(2-K5)*4 

J  6  =  ( K  5  - 1  1  •tNM2 

00  44  .1  =  1  .NT  .N5 
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U1  =  0. 

J1=C J-l )*NM4+J5 
DO  47  1=1, NM? 

J2=Jl+I 

IJ1=U1+VR(  J?)*F3(  I  ) 

47  CONTINUE 
J8= J7+J 
E( J8)=C1*U1 
S2=CABS(Et J8) ) 

51  fit J8)=S2*S? 

52  =  SI  G (  J8 )  >:<SCL  (  K4  ) 

E9{ J)=b.+SN( J)ms? 

C'/W  I  t  ~C  «  t*  »  i  »  >, 

u  tut  O/'-sX 

44  CONTINUE 

DO  86  J  =  1 ,  N  T  » N S 

CALL  SYMBOL  ( E9  ( J  )  ,  b  1 0  (  J  )  ,  *  07  t  K2 , 0.  , -1  ) 

86  CONTINUE 

00  89  J  =NT 4 , NT 3, NS 
J 1=NT- J+  1 
Sl=10„-E9{ Jl) 

CALL  SYMBOL! SI, E10( Jl ) , .07,K2,0. ,~1 > 

89  CONTINUE 

CALL  L I N E ( E 5 ( K 3 ) ,E6 < K3 ) ,NT2, l , 0, 0 ) 

66  CONTINUE 

68  CALL  LINEUP, YP, 2, 1,0,0) 

DO  77  J  =  1 ,7 
Sl=9-J 

CALL  SYMBOL (SI, 5. , .14,13,0. ,-l ) 

77  CONTINUE 

CALL  LINE(YP,XP,2,1,0,0) 

DO  73  J= 1 , 7 
Sl=9-J 

CALL  SYMB0L(5.,Si,el4,13,90.,-l) 

78  CONTINUE 

CALL  PLOT ( 7 . , 0. ,-3 ) 

WRITE (3 , 122 )  K 

122  FORMAT ( *  0* , 1 3 ,  *  MODE  SCATTERED  FIELD  AND  SCATTERING  CROSS  SECTION/ 
1WAVELENGTH  SQUARED1) 

Sl=(2-K4m80« 

WR I TE ( 3 , 1 29 )  SI 

129  FORMAT! (  INCIDENCE  FROM  0=»,F4„0,«,  0=0') 

WR I TE ( 3 , 1 20 ) 

120  FORMAT ( 5  + 1 , 15X , 1 - 1 , 7X , 1 / ' ) 

WRITE (3,39) 

WRI TE(3? 117 ) 

DO  126  J=i,NT,NS 
J 1 = J+NT 

WRITE (3. 38)  TH(J) tE( J) tSIGf Ji*E( Jl)fS!G( Jl) 

126  CONTINUE 
65  CONTINUE 
29  CONTINUE 

CALL  PLOT ( 6 . , 0 . , -3 ) 

STOP 

END 
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VII,  GAIN'  PATTERNS 


Program  it 6  which  calculates  and  plots  the  gain  pattern  for  radiation 
from  axially  symmetric  excitation  on  a  surface  of  revolution  accepts  punched 
card  data  according  to 


10 

11 

15 

50 

57 


READ (1,10)  NT ,  NT,  NS,  JM,  NC,  NV,  BK,  SCL 
FORMAT (613,  2E14.7) 

READ(1,11) (AMD(I) ,  1=1,  JM) 

FORMAT  (5E14.7) 

READ (1,15) (RH(I) ,  1=1,  NP) 

READ (1,15) (ZH(I) ,  I  *  1,  NP) 

FORMAT  (10F8.4) 

READ(1,50)(L(I),I  =  1,  NC) 

FORMAT  (2013) 

READ(1,57) (V(I) ,  1=1,  NV) 

FORMAT  (7E11.4) 

READ (1,50) (LV(I) ,  1=1,  NV) 


The  variables  NP,  NT,  NS,  JM,  BK,  AMD,  RH,  ZH,  and  L  are  the  same  as  in 
program  it  5  except  that  L  now  applies  only  to  the  directed  eigencurrents . 
The  output  of  program  H2  characterizes  each  eigencurrent  by  NM2  numbers. 

For  the  axialxy  symmetric  mode,  either  the  first  NM2/2  numbers  (correspond¬ 
ing  to  J^)  or  the  last  NM2/2  numbers  (corresponding  to  J^)  are  supposed  to 
be  zero.  Thus  L(I)  indicates  that  the  L(I) ^  eigencurrent  (only  ut  directed 

i.t.> 

eigencurrents  being  considered)  to  come  out  of  program  it 2  will  be  the  Icn 

current  to  be  added  to  the  modal  expansion.  There  are  NV  _<  NM  axially 

rh 

symmetric  slots  on  the  body  of  revolution.  The  voltage  across  the  I  slot 
is  V(I) .  Note  that  V(I)  is  complex.  The  slot  is  located  at  the 
(2*LV(1)+1) data  point  (p  =  RH,  z  =  ZH)  on  the  generating  curve  C.  Accord¬ 
ingly,  the  IL^  slot  occurs  at  the  peak  of  the  LV(I)t^1  triangle  function 
°^LV(i)‘  ‘^e  -unc"or*  appears  in  (2-42)  of  [!}.  For  the  plots,  one  inch 

corresponds  to  a  gain  of  l./SCL. 


49 


The  impedance  matrix  and  eigencurrents  are  read  from  the  third  and 
fourth  records  of  direct  access  data  set  6  according  to 

REWIND  6 
READ  (6) 

READ  (6) 

READ  (6)  (Z(I),  I  =  1,  NZ) 

READ  (6)  (FI (I),  1=1,  NZ1) 

tt 

Only  half  of  the  impedance  matrix  is  read  because  only  [2^  ]  of  (2-46)  of  [1] 
is  needed.  Program  //I  has  stored  the  impedance  matrix  on  the  third  record  of 
data  set  6  columnwise  in  the  block  diagonal  form  given  by 

[z“]  [0] 

fz  ]  «  (22) 

o 

to]  [zjh 

L  O  J 

The  subroutine  PLANE  is  similar  to  the  one  compiled  with  program  #5. 

^  t  0 

At  present,  only  IT  of  (5)  is  needed  because  a  slot  voltage,  corresponding 
to  a  ufc  directed  axially  symmetric  aperture  electric  field,  induces  only  a  u^ 
directed  electric  current  which  radiates  only  a  u^  directed  far  field* 

DO  loop  60  suppresses  the  lower  left  zero  submatrix  on  the  right  hand 

side  of  (22).  The  net  result  of  DO  loops  71  and  73  is  to  arrange  the 

eigencurrents  (FI)  and  their  eigenvalues  (AMD)  in  the  order  dictated  by  L 

and  to  suppress  intervening  zeros  from  FI.  DO  loop  62  obtains  the  matrix 

[X]  appearing  in  (2-15)  of  [1]  by  taking  the  imaginary  part  of  the  average  of 

corresponding  off  diagonal  elements  of  the  impedance  matrix.  Statement  79 

inverts  the  impedance  matrix  to  obtain  the  admittance  matrix.  DO  Icop  82 

obtains  the  column  matrix  [Iq]  associated  with  the  electric  current  by  pre- 

tnultiplying  the  excitation  column  matrix  by  the  admittance  matrix.  The 

elements  of  the  excitation  column  matrix  of  (27)  of  [3]  are  actually  2n  times 

the  slot  voltages,  but  the  factor  2tt  is  inconsequential  as  far  as  that  gain 

£ 

is  concerned.  DO  loop  82  stores  [I  ]  in  E3.  DO  loops  84  and  97  compute  the 

1 6  t  ^ 

radiation  field  E/~  R  [I  ]  and  gain  and  store  them  in  E  and  G.  The  radi- 
8  o  o  *  8 

ation  field  and  gain  are  normalized  so  that 


f'T 
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I  G,  sin  t*  d0  =  2  (23) 

j 

0 


The  phase  factor  -je  suppressed  from  E^.  DO  loop  99  prepares  hori¬ 

zontal  and  vertical  components  E5  and  E6  suitable  for  plotting  the  gain. 

DO  loop  48  multiplies  the  eigencurrents  by  p  to  retrieve  the  column  matrix 
[L^]  associated  with  each  eigencurrent .  The  index  J  of  DO  loop  48  indicates 

Lite  eigencurrent.  DO  loop  21  stores  the  matrix  [I] - - - 

[X] [X] [I]  (j  +  1/A) 

appearing  in  (20)  in  T3. 

The  index  K  of  DO  loop  29  indicates  that  inner  DO  loop  31  will  add  the 
contribution  from  the  K^1  eigencurrent  to  the  admittance  matrix  (20).  DO 

t  A  t  G 

loop  44  premultiplies  [1^]  stored  in  E3  by  Rq  stored  in  VR  to  obtain  the 
radiation  field  and  the  gain.  The  radiation  field  EA  and  gain  GA)  normalized 
according  Co  (23),  are  stored  in  E  and  G.  In  order  not  to  mask  a  possible 
discrepancy  in  the  amplitude  of  the  approximate  pattern  obtained  by  super¬ 
imposing  a  few  eigenfields,  the  present  normalization  uses  the  value  of  the 
pattern  integral  (23)  previously  computed  from  the  admittance  matrix  obtained 
by  inverting  the  impedance  matrix. 

Minimum  allocations  are  given  by 

COMPLEX  Z(NM2*NM) ,  Y (NM2*NM) ,  VR(NT*NM) ,  E3(NM) , 

T3(NM*NM) ,  E(NT) ,  V(NV) 

DIMENSION  AMD ( JM)  ,  RH(NP) ,  ZH(NP) ,  FI(NM2*JM), 

DH(NP-l),  TH(NT) ,  T2(NM*NM) ,  X(NM*NM)  , 

G(NT) ,  L(NC) ,  LV(NV),  SN(NT),  CS(NT),  GX(NT) , 

GY (NT),  E5(NT*2),  E6(NT*2) 

COMMON  RS(NP-l),  ZS(NP-l),  SV(NP-l) ,  CV(NP-l) ,  T(NM*4) 

DIMENSION  3J(3*(NP-1)) 

DIMENSION  LR(NM) 

where 

NM2  =  NP  -  3 
NM  »  NM2/2 

Here,  NT  is  its  value  after  execution  of  statement  90  in  the  main  program- 
NY- 1 e  that  BJ  appears  in  PLANE  and  LR  in  LINEQ. 
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Listing  of  Program  #6 

//  (  0034,  EE,  3,?,  ,6) ,  'MAUTZ  ,  JOE  '  ,MSGI.EVF.L=  1 

//  EXEC  FORTGCLG » P ARM  .  Fl)RT  =  1  MAP  ' 

//FORT.SYSIN  00  * 

SURR0U1 I  OF  PLANE (VVR  ,THR  ,NT  ) 

COMPLEX  VVR  (  1  )  , Ab,A6,U 

COMMON  U ,  R  S  (  4  0 ) ,7S(40) ,SV(40) ,CV(40) ,BK,NP,T(HO) 
DIMENSION  R  J ( 126) , THR ( 1 ) , FK ( 20 ) 

NN  =  0 
KG=NP-1 
NM=KG/2-l 
M2=NN+2 

A5  =  2  .*3 ,141  59  3*U*#  ( NN+ 1 ) 

F  K  (  1 )  =  1  , 
no  153  J=l,M? 

J1=J+1 

FK ( J 1 ) =FK { J  )  #  J 
153  CONTINUE 

DO  155  L=1,NT 
L  1  =  ( 1.-1  )  +NM 
CS=C0S(THR«L) ) 

SN=SIN(THR(D  ! 

BCS=BK#CS 
no  302  J  =  1 , KG 
X-RS( J)*BK*SN 
J1  =  J 
1 1=NN 

IF  (ID  303,304,303 

304  11=11+1 
J1=J1+KG 

303  DO  305  JJ=11,M2 
I F ( X-l . E-5 )  1,1,2 

1  IF(JJ-l)  3,3,4 

3  RJ ( J 1 >  =  1  . 

GO  TO  306 

4  R.J  ( Jl ) =0. 

GO  TO  306 

2  RH=X/2. 

RH.2  =  RH*RH 
RH3=RH**\ JJ-l) 

RJ(J1)=RH3/FK(JJ) 

SS=RJ( Jl) 

8  SST=SS*l.F-7 
DO  155  K= 1 , 20 
SS=-SS*RH2/K/(K+JJ-1 ) 

R J  (  Jl )  =  RJ  ( J 1 )  i-SS 
I F ( ABS ( SS ) -SST )  306,306,155 
155  CONTINUE 
STOP  155 

306  J1=J1+KG 

305  CONTINUE 
302  CONTINUE 

IF(NN)  307,308,307 

308  DU  309  J=l ,K0 
J1=J+2*KG 

BJ  ( J )=-RJ ( Jl ) 

309  CONTINUE 

307  DO  300  J=1,NM 
J1=J+L1 

VVR ( Jl )=0. 

DO  301  1=1,4 


11=2- ( J-l 5+1 
Ift=9>-(  j-1 5 1 ! 

I  2= i 1+KG 
; 3= i 7+kg 

A*  =  cc»isf  7SI  1 1 »' HCS1+U*SIN(ZS(  1 1  )*HCS)  )<A5 
h  j2=  ( b.i (  m-oji  1 1 )  >*.s 

VVR(  J1)=VVE(,U  )+A<s*(CS<-SV<  1 1 ) *BJ2+SN$CV (  I2)*U)*T(14) 

HOI  CONTINUE 
300  CtHfl'HUe 
ISA  CONTINUE 
Pi-TURN 
K  I'M 

SUOfiljilT I Nt-  I.INE  90.1,0 
COMPLEX  Cll)  .STliR  , STO , ST  ,  S 
OMENS  I  ON  L«fH8) 

•Ki  20  1  =  1, LL 
LR<I)=I 
?n  CONTINUE' 

Mj  =0 

on  io  w=i,ll 

K=M 

on  2  i=m,li 

K  1  =  M )  +  I 
K?=M1+K 

I E { C AOS ( C ( K 1 5 )-CAHS(C(K2) ) )  2, 2,  ft 
8  K  =  I 

7  CONTINUE 
LS=LR(M) 

IR{M)=LR(K) 

LR(K)=LS 
K2=M1+K 
S  TOR  =C ( K2 ) 

Jl=n 

00  7  J  = 1 ,  L  L 
K1=J1+K 
K?=J1+M 
sm=c(Ki ) 

C (K1 )=C(K?) 
c(K2)=srn/sTnR 
J1=J1+LL 
7  CONTINUE 
KlsMl+M 
C(K1)=1./STDR 
01)  11  I  =  1 » LL 
IE(I-M)  12.11,12 
1?  K 1 sM 1 + I 
ST  =  C ( K  1 ) 
f.(Kl  )  =0 . 

J1=0 

on  10  J=1 , LL 
K  1  =  J  1  +  I 

!;?=ji+k 

C(K1)=C(K1)-C(K2)*ST 

J 1  =«J  1+LL 

10  CONTINUE 

11  CONTINUE 
M  =M1*LL 

1«  CONTINUE 
.11=0 

Ol>  9  J=l, LL 
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IF(J-LR(J)J  14, P.14 

14  LRJ=LR(J) 

J2  =  (LRJ-1 )  '-LL 
21  Dl)  13  1  =  1,  LL 
K2=J?+ J 
K1=J1+1 
S=C ( K? ) 

C(K2)=C(K1 ) 

C(K1)=S 
13  CONTI NOE 

LRi J)=LR(LRJ) 

L«(LPJ)=LR.I 
IRJ-LR(J))  14,8,14 
k  Jl=Jl+LL 
9  CONTINUE 
RETURN 
END 

COMPLEX  11,111  ,Z  (  7??)  ,YU4*4  J  ,VR  ( 1387)  ,E3(  19)  ,7 3( 3M  )  »E ( 73) , V(  14>) 
DIMENSION  AM|i(38),RHI41),ZH(41),FM  1444  )  ,  OH{  40  J  ,  TH(  73  )  ,  T2  (  361  ) 
f)  1 PENS  1 ON  X( 341)  ,0(73) , L ( 14  1 , LV ( 19 ) , AREA ( 400) ,SN(73) ,CS(73) 
DIMENSION  GX{ 7^)  ,GY (73 ) ,XP(2) , YP(2 ) ,E5( 146) ,  EM 146) 

EMU  I  VALENCE (T2(l 1 , X (1 ) ) , (2 ( 1 1 , Y U ) ) 

COMMON  U,RS(40) ,75(40) ,SV(40) ,CV(4H) ,BK,NP,T(HO) 

CALL  PL0T5( AREA, 400) 

READ! 1, 10)  NP *NT,NS » JM,  NC , MV , BK, SCL 

10  FORMAT ( 61 3, 2E 14. 7) 

RE AD( 1 » 1 1 ) ( AMD ( I ) , I = 1 • JM ) 

11  FORMAT ( 5E 14 .7 ) 

RFAiM  1,15) (RH( I  1,1  =  1, MP  > 

READt 1,15) (ZH( 1 ) .!=) ,NP) 

15  FORMAT (10F8. 4) 

REALM  1,50) (L(l) . 1  =  1, NC» 

50  FORMAT (201 3) 

READ ( l , 57 ) ( V i () , I = 1 ,NV ) 

57  FORMAT (7E 11.4) 

R6AD(l,50)(LV(n,r=l,NV) 

WRITE (3, 37) 

37  FORMAT ( • 0  NP  NT  NS  JM  NC  NV ' ,6X , * BK » , 12X , ' SCL » i 
W'R  I  TE  (  3, 3ft  )  NP  , NT  , M S ,  JM , NC , N V ,  HK ,  SC L 
3H  FORMAT ( 1X,6I3,2E14.7) 

-WRITE  (3,5)  ( AMIM  1 ) ,  1  =1 ,  JM) 

5  FORMAT ( •0AMD’/(1X,5F14,7) ) 

VRITE(3,39)  (RH(  M,I=1,NP; 

39  cORMAT( *  ORH  * / ( IX , iOFR .4 ) ) 

WRITE  (3, 9  )  UHl  I)  ,  I  =  1,NP) 

9  FORMAT  (  1  07.H  1  /  (  IX  ,  lOFft  .4  )  J 
WR I TE ( 3  *  8  > ( L ( I ) , 1  =  1, NC) 

8  FORMAT ( :0L' /( 1X.20I 3) ) 

WRITF(3,7)(V(l),i=l,NV) 

7  FORMAT l  ‘  OV 1  / ( LX » 7E 1 i  .4  } ) 
ut<  j  tp  {  3 .  a  }  (  LV!  ! ) :  1  =  1  ,NV ) 

A  FORMAT  (  '  OLV  1  /  (  IX ",  201 3 )  ) 

PI=3. 141593 
0=10. ,1.) 

IF( (RH(I)-RH(NP) ) .NE.O. .OR. (ZH( l )-ZH(NP> I .NE.O. )  G.)  TO  90 
RHIMP+l )=RH( 2 ) 

ZH(NP+1 )=7K( 2) 

RH ( NP+2 ) =RH( 3 ) 

ZH(NP+2) =70(3) 

NP=NP+2 


90  »iM2=NP-3 
NM=NM2/2 
Nr=NM?*MM 
N7 1=NM?*JM 
REWIND  ft 
READ(ft) 

RFAf)( 6) 

READ (6) ( Z  t I ) ,  1  ~1 ,N7 ) 

READ(ft) (FI ( I ) , I = 1 ,NZ 1 } 

NT  ?=NT  *2 
NT3=MT-N<; 

NT4=NS+1 
DU  40  I =?  ,NP 
12=1-1 

P.R1=RH{  I  )  — RH (  I?) 

RR2=ZHll  WHU?) 

DM ( 1 2 ) =SORT ( RR 1 *RR 1 fRH2#RR2) 

RSU2)  =  .5*(R“(  I  )+Rh;  I?) ) 
ZSU2)=.5*(ZH(  D+ZHU2J » 
SV(I?)=RR1/DH{ I?) 

CV( I2)=RR?/DHl I?) 

40  CONTINUE 

I  HI  41  J=  1 « MM 
J2=2* ( J-l ) + 1 
J3=J2+1 
J4=J3+1 
J5= J4+ l 
J6=4*{  J-1JM 
J7=J6+1 
J8= J7+1 
-J9= Jfl+1 

DEL1=DH  ( .1? )  +DH I J3 ) 

DEL2=DH{ J4 1 +DH  ( J5) 

T  ( J6 )  =DK(  J?)  *DH(  .12 )  /  2»/DELl 
T(J7)=DH{J3)-(DH(J2) +DH (J3)/2« ) /DEL l 
T  (.18 )  =DH ( .14 )  - ( DHI .15 )  +DH l  J4 )  /2 .  )  /DEL2 
T( J9)=DH( J5)*DH( J51/2./0EL2 

41  CONTINUE 
DEL=PI /(NT-1 } 

Df)  43  J  =  1 ,  NT 
TH( J )= ( J-l J -DEL 

SN( J1=SCL*5IN(TH( J!) 

CS( J)=SCL*COS(TH( J) ) 

43  CONTiNUF 

CALL  PLANE ( VR , TK,NT ) 

XP ( 1 ) =2  . 

XP( 2) "ft. 

YP(1 5=5. 

YP{ 2)=5. 

DEL 1=180. /PJ 
Oil  64  J= 1 ,NT 
TH(  j)=TH(.|)*dEL1 
M  CONTINUE 
J1=0 

DO  60  .1=1,  NM 
.l?=(.J-l)*NM2 
DO  61  1=1, NM 
Jl=Jl+l 
,13= J2+I 
Z(J1)=Z(.I3J 
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61  CONTINUE 

h 0  CDNTINUr- 

J1=0 

DO  71  J=l,NC 
J?.*CL(J)-1)#NW? 

DO  7?  1=1  ,NM 
Jl'-Jl  +  l 
J3= J2+ I 
121 J1 l-FI! J3> 

72  CONTINUE 
J2=L(J) 

E5(J)=AM0(,I2) 

71  CONTINUE 
JT  =0 

V  3  J=1,NC 
on  74  1=1 ,NM 
Jl=Jl+l 
FI(J1)=T2U1) 

74  CONTINUE 

AMD(  .1  I="E5  ( .1) 

73  CONTINUE 

00  62  J=1,NM 
J2=( J-1J0MM 
on  63  i=i,. i 
Jl=J2+ I 
J3=( 1-1 )*NM<J 

Xt Jl)s.5*AiMAG(Z(Jl)+ZCJ3>) 

XI J3)=XIJ1) 

63  CONTINUE 

62  CONTINUE 

79  CALL  LINEOINH.Z  I 
00  82  1=1, NM 
E3( I ) =0. 

DO  6b  K  =1 ,NV 
K1=I+(LV(K)-1} "NM 
E3  I  I  )  =E3  I  I  )  +7  t  K 1  )  <=V  (  K  ) 

65  CONTINUE 
82  CONTINUE 

S2=0. 

J?=0 

DO  84  J=1,NT 
U1=0. 

DO  88  1=1, MM 
J2=J2+1 

U1=U1+VR ( J? )  **E3  ( I ) 

88  CONTINUE 
El Ji=Ul 

si=cabs;ui) 

G( J)=S1*S1 
52»S2+G(J)«SNU) 

84  CONTINUE 

S5=2.*SCL/S2/DFL 
S6=S0RT(S5J 
DO  97  J=1,NT 
EIJ)=E(J.)*S6 
G  ( J )  =G  (J )  *S5 
97  CONTINUE 
WR1TE(3,66) 

66  format i *  or ad i at  ion  field  ano  gain  by  matrix  inversion*) 

UP.  ITE 1 3 ,67  ) 


67  FORMAT  (  1 0  0  REAUEO)  I  MAG  (ECU  GAINO'  ) 

WR I TE ( ?  ,  68 ) 

68  FORMAT ( ' +  -• , I1X, f 1 1 X , »-* ,10X, •-» ) 

on  69  vl  =  l.'!T,N9 

WR  1 1 F  {  3 »  7  0 )  TH{J),F(J),G(JJ 
70  FORMAT ( 1 X  ,  86 . 1 , 3F  1 2 .4  ) 

69  CONTINUE 

00  99  JsltNT 
Sl=SN( J)»G( J) 

J9=NT2- J+ 1 

E5( J)=5.+51 

E5 I J9 ) =b.-S 1 

E6( J)=5.+CS( J)*G( J) 

E6( J9)=E6( JJ 
99  CONTINUE 
Jl=0 

nil  48  .J=1,NC 
00  18  I=1,NM 
J1=J1+1 
J4=2*I+1 

FI  ( jn=FHJl)*RH{J4) 

18  CONTINUE 
48  CONTINUE 

00  21  J=1,NC 
Jl={  J-1)*NM 
S1=0. 

1)0  22  1=1, NM 
S2=0. 

J4=( I-1)*NM 
1)0  23  K=1»NM 
J2-J4+K 
J3=J l+K 

S2=S2+X( J2)*FI ( J3) 

23  CONTINUE 
J3=J1+I 

Sl=Si+S2*FI( J35 
22  CONTINUE 

111  =  1. /Sl/IU+l./AMDIJ)) 

DO  27  1=1, MM 
J2=J1+I 

T3(  J2)=FI  ( J2)!)[U1 
27  CONTINUE 
21  CONTINUE 
NZ">=NM*NM 
DO  32  J= 1 »NZ 2 
Y  (  J  )  =0  . 

32  CONTINUE 

nri  29  K  =  1  ,NC 
J3=  (K-l )  =>NM 
J1=0 

00  31  I=i,NM 
J4= J3+ 1 
E3 ( ! ) -0. 

DO  75  J=1,NV 
J2= J3+LV { J ) 

J1=J1+1 

Y(Jl)=Y(Jl)+T3(J/*)<tFI  (J2) 

F3U)=E3I  n+Y(Jl)*VIJ) 

75  CONTINUE 
31  CONTINUE 
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1)0  44  Js  1. NT, NS 
l'1-O. 

ji=(  j-h^nm 

DO  47  1=1,  NK 
J2= J1+ I 

U1=U1+VRM?)*E3(  I  ) 

47  CONTINUE 
ElJ)=Ul*SA 
S1=CABS(E( J) ) 

G< J)=S1*S1 

GX( JJ=5.+SM( J)*G( J) 

GY ( J } -5  *  +CS ( J ) *G ( J ) 

44  CONTINUE 

WR I  T  E  ( 3 , 7  A )  K 

7 A  FORMAT! 'O' , 13, '  MODE  RADIATION  FIELD  AND  GAIN*  j 
KR I TB ( 3 , A7 ) 

WR I  TE  (  3 » A8 ) 

DO  80  J=1,NT,NS 

WRITE  I  3, 70)  TH( J) ,E(J) ,G( J) 

80  CONTINUE 

CALL  LINEIXP, YP, ?, 1,0,0) 

DO  77  J=1 ,7 
Sl=9-J 

CALL  SYMBOL ( S 1 ,5  . « . 14, 13,0. ,-1) 

77  CONTINUE 

CALL  L I NE ( YP , XP , 2 , 1 , 0  *  0 ) 

DO  78  J  =  1 , 7 
Sl=9— J 

CALL  5YMB0L(5.,S1,.I4.13,90.,-I) 

78  CONTINUE 

CALL  LINE(E5,E6, NT?, 1,0,0) 

DO  8A  .1=1, NT, NS 

CALL  SYMBGL(GXIJ) ,GY(J) ,.07,4,0., -1) 

8 A  CONTINUE 

DO  89  J=NT4,NT3,NS 

J1=NT-J+1 

Sl= I0.-GXI J1 ) 

CALL  SYMBOL  (SI, GY  (  Jl)  ,  .07  ,4 ,0.  ,->l ) 

89  CONTINUE 

CALL  PLOT ( 7. ,0. ,-3 ) 

29  CONTINUE 

CALL  PLOT (6. ,0., -3) 

STOP 

END 

/* 

//GO.FTOAFOOl  DD  0SNAME=EE0034 .REV1 ,DI SP=QLD,UN IT=2314 ,  X 

//  VQLUME=SER=SU0004 ,DCB= ( RECFM=V »6LKSIZ£=1 800, LRECL =1796) 

//GO. SYS IN  DD  * 
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